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TECHNICAL MEMORANDUM 


EQUATIONS OF MOTION OF A SPACE STATION WITH EMPHASIS 
ON THE EFFECT OF THE GRAVITY GRADIENT 

SECTION 1 - MODEL DESCRIPTION 


The class of space stations treated in this report does not include space sta- 
tions with the most general configuration conceivable. The station model dealt with 
herein is comprised of a rigid central body supporting several flexible appendages, 
smaller rigid bodies, and point masses. The smaller rigid bodies include swivel 
engines and rotors while the point masses may be part of the mechanical analog of a 
consumable liquid aboard the carrier or may represent small bodies having a pre- 
scribed motion relative to the carrier (trim masses, for example, installed for the 
purpose of attitude control) . Each flexible appendage admits 3-D bending and has 
at most two "rigid body" rotational degrees of freedom relative to the central body. 
The flexible appendages are not interconnected. The disturbing effects of crew 
motion on vehicle attitude will be ignored. 

This particular choice of space station model was strongly influenced by the 
author's brief acquaintance with the SEPS vehicle of the early 1970's to which he 
applied certain of his previous developments [1] in describing its motion. 


SECTION 2 - COORDINATE SYSTEMS 


The formulations of this paper required the introduction of the following 

rectangular coordinate systems. 

Vn z n — The N -frame is a Newtonian frame whose existence is postulated (other- 
wise there would be no Newtonian mechanics). It is stipulated here that 
the N- frame is oriented as the S- frame defined below (this being merely 
for convenience, actually). The N-frame will enjoy brief recognition, its 
introduction being made solely for the purpose of preventing failure to 
include certain terms in the equations of motion (in particular, the vector 
translational equation). See Figure 1. 

XgYgZg — The S- frame has origin at the geometric center of the reference ellipsoid 
with the Xg axis directed through the mean vernal equinox of a specified 

epoch, the Zg axis directed as the Earth's mean spin vector, and the Yg 

axis (lying in the mean equatorial plane) so directed as to make the sys- 

tem right-handed. The unit vectors Ig , J g , and Kg which span the 

S-space are directed as the Xg, Yg, and Zg axes, respectively. The 

S-frame is coordinate system No. 4 of Reference 2. 

XgYgZg — The origin of the "Earth fixed" or E-frame coincides with that of the 

S-frame and the Zg axis is directed as the Zg axis. The XgYg plane is 

the mean equatorial plane, the Xg axis being directed through the prime 



meridian (for all t >. t ) and the Yg axis directed so that XgYgZg is 
right-handed. 

The B-frame (the structural axes system), also right-handed, is at rest 
relative to the rigid central carrier, and its origin and orientation may be 
arbitrarily chosen relative to the central body. Notice that this definition 
does not require that the origin of the B-frame be embedded in the cen- 
tral body, though such may be the case if deemed more convenient. (In 
the early Saturn vehicles, for example, the structural axes had origin 

► 

100 in. behind the engine gimbal plane.) The unit vectors i, j, and k 

spanning the B- space are directed as the x, y, and z axes, respectively. 
See Figure 1. 

The "body" axes system (herein designated the B-frame) has origin at 
the instantaneous center of mass (CM) of the entire vehicle system and 
is oriented always as xyz (the B-frame). The origin of xyz has position 
vector fgjyj relative to xyz. As a consequence of the redistribution of 
of mass due to bending, motion of internal parts, etc., the vector r^ 
will vary in both direction and magnitude. The unit vectors i, j, and k 
also span the B -space. See Figure 1. 

The i- frame, a right-handed system, has origin at the "idealized" point 
of attachment of the i in flexible appendage, the x^ axis serving as the 
axis of rotation about which flexible appendage i can rotate (through the 
angle e.) relative to the rigid central body. Bending displacements of 
points belonging to flexible appendage i are referred to axes x.y.z., that is, 
by the symbol A 1 = A (r^, t) is meant the displacement due to bending of 
the point which had position vector r. (referred to x.y.z.) before deforma- 
tion. The unit vectors i., jj, and kj span the i-space (i = 1, ..., NA), 
these being directed as the x^, y., and z^ axes, respectively. It should 
be remarked here that the attachment point, being a point of the central 
body, does not itself undergo a "bending" displacement by virtue of the 
assumed rigidity of the central body. See Figure 1. 

These axes have origin at the instantaneous CM of flexible appendage i 
whatever its configuration, deformed or undeformed. They are always 
oriented as x,.y.z. (i = 1,... NA). 



x Ri y Ri z Ri — Ri"f rame has origin at the CM of rotor i (i = 1, ... NR), the x R . 

axis being coincident with the axis of rotation and directed in accordance 

'-V~ 

with the right-hand rule. The Ri-frame is fixed relative to the B -frame. 

It too is a right-handed system as are all the other systems of this report, 

— ► — > — > 

Spanning the Ri- space are the unit vectors i R j, j R j, and k Ci . . 


and k T 


X Ri y Ri z Ri ~ The is fixed relative to rotor i and coincides with the 

Ri-frame at the instant the rotor begins to rotate and at such times at 

which the rotor completes a revolution. (Introduction of x D !y_ !z D ! was 

necessary because of the method of formulating the equations of motion. 

The angle, through which rotates relative to 

will not appear in the final results though its first and second time 

derivatives will.) The (Ri)' space is spanned by the unit vectors h,.’, 

Kl 

W and E Ri'- 

x Gi y Gi z Gi ~~ T ^ e &™^ a ^ frame axes (Gi-frame) of the i** 1 single DOF (SDOF) control 
moment gyro (CMG) have origin at the CM of the gimbal frame, the x^^ 
axis being coincident with the gimbal axis of rotation and the z^,. axis 
directed as the spin vector of the gyro element. The x Q ^ and z^. axes 
are directed in accordance with the right-hand rule. The direction of 


the x_. axis is invariant in the B- and B- frames. 
Gl - 

jg., and k^,. span the Gi- space. 


The unit vectors i. 


x .y 

ffr 


gh gi gi 


The gyro element axes of the r n SDOF CMG (the gi-frame) have the same 
origin as ^^G^Gi’ ^ being assumed here that the gimbal and gyro 
element have a common CM. The z ^ axis is directed as the gyro element 
spin vector. The gi-frame is fixed relative to the gyro element and coin- 
cides with the Gi-frame at each instant it completes a revolution. The 
— ^ — > — ► 

unit vectors i ., j^., and span the gi-space. (As with the x R !y R ! z R ! 
frame it was the method of formulating the equations of motion that 
required introducing the x .y .z . frame. The angle v •, through which 

&DD « & 


the gi-frame rotates at the constant rate v? . = w . relative to the Gi- 

& gi gi 

frame, will not appear in the final results, but to . will be present in 
several terms.) 


Associated with a two degree of freedom (2DOF) CMG are the following four 
coordinate systems and corresponding vector bases. All are right-handed and have 
the same origin, namely, the CM of the outer gimbal-inner gimbal-gyro element com- 
bination (it having been assumed that those three components of the CMG have a 
common CM.) 


3 



X GB^GB Z GB — The base frame (GB-frame) is fixed relative to the B-frame and 

coincides with the OG -frame when the outer gimbal angle 6 QG = 0. The 
GB- space is spanned by the unit vectors i^ B , T GB > and k. 


GB * 


x OG^OG z OG _ The ou ^- er gimbal space (OG- space) is spanned by the unit vectors 

Iqg ’ 1qq> anc * ^OG' The z OG axas (directed as kg G ) is the outer gimbal 
rotation axis and its direction is invariant in the GB- and B- frames. The 
OG-frame is oriented as the IG-frame when the inner gimbal angle 6 j G = 0. 

— ► 

X IG^IG Z IG _ The ainer gimbal space (IG-space) is spanned by the unit vectors ij G , 
■*IG’ an< * ^IG ’ The x jQ ax i s (directed as i^ G ) is the inner gimbal rotation 
axis. The gyro element spin vector is directed as . 

(The unit vector triads i j k , E, - GB , OG, IG , above are those of Reference 3 in 
different notation.) s 


x y o* ~ The gyro element space (g-space) is spanned by the unit vectors i , j , 
and kg. The y^ axis, directed as both j and j IG , is the gyro element 
spin axis. The g- frame (fixed relative to the gyro element) coincides 
with the IG-frame initially (that is, at the time designated t ) and at each 
instant it completes a revolution. 


(The seemingly absurd introduction of the axes x ,y z , never used in any of the 

g g g 

literature browsed by the author, is necessary to the formulations of this paper, 

although the angle <p through which the g-frame rotates at the constant rate <i > , = w 

g g g 

relative to the IG-frame will not appear in the final results. The spin rate <*> , how- 

g 

ever, will appear in certain terms of the moment equations, the outer gimbal equation, 
and the inner gimbal equation.) 

The following is associated with the i** 1 swivel engine: 


x Ei y Ei z Ei — The El-frame has origin at the CM of swivel engine i, is fixed relative 
to swivel engine i, and if sensibly oriented has one of its axes directed 
as the engine thrust vector. The Ei- space is spanned by the unit 
vectors i^., j^., and kg. directed as Xg^, y^j. and Zg^, respectively. 
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SECTION 3 - TRANSFORMATIONS AND SUBSIDIARY RELATIONS 


The vector bases defined in Section 2 give rise to the rotation matrices* defined 

below with important attendant relations. 

T — The upper case letter T (when not used as a superscript to indicate the trans- 
pose of the matrix or column vector to which it is attached) will denote the 
rotation matrix defining the transformation T(S-*-B), the literal translation of 
the symbol T(S->-B) being "transformation of the resolution of a vector on the 
S- vector basis to the resolution of that vector on the B- vector basis," or, more 
succinctly, "transformation from the S resolution to the B resolution." From 
the definitions of the B- and B-frames, it should be evident that T(S->-B) and 
T(S-*-B) are completely equivalent. (Insofar as consistency in notation is con- 
cerned, some will be eager to point out, in view of certain definitions to follow, 
that better notation might have been realized by attaching the subscripts SB to 
T . Though the author will be the first to agree that attaching the subscripts 
SB to T would, for the sake of consistency, be better, his frequent use of the 
letter T without subscripts in both years past and recently has left him almost 
impervious to change.) 

The matrix T cannot be completely specified until an "Euler sequence" has been 
decided, that is, until one prescribes the sequence of rotations which the S- 
frame would have to undergo to assume the orientation of the B-frame (and 
hence of the B-frame). For a two, three, one sequence through the angles 
y P n , </>. • , and <p , respectively, T is given by 


T [ *r ] (l) [ V(3) [ *P ] (2) 


where 


[ *P ] (2) 


cos *> p 
0 

_ sin <p 


0 -sin <p 

r 

1 . 0 

0 COS y?p - 


(3-0.1) 


*A11 rotation matrices herein are "proper" rotation matrices. 
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[ V(3) 


COS ip 


-sin y? 


sin 


cos 


0 

0 

1. 


(3-0.2) 


[ *r ] (l) 


1. 

0 

0 


cos 


■ sin < p 


sin v 


cos J 


(3.03) 


Following the method laid down in References 4 and 5, among others, it is shown in 
Reference 1 (and no doubt elsewhere) that corresponding to the Euler sequence above 
the angular velocity, uTg = [u>^, wg, u^] , of the B-frame (or B-frame) relative to 

the S-frame is given in terms of the Euler angles and their first time derivatives by 



- U1 - 


sin <p 

y 

0 

l.“ 


1 

•S* 

I 

111 

w 2 

— 

COS i fi COS < fi 

y r 

sin 

0 


fi 

y 


- w 3- 


_ -cos sin fi p 

cos fi r 

0 _ 


i 


(3.04) 


The vector wg as given above has the B (and B) resolution. It is also shown in 
Reference 1 (and in literature not cited) that 


E T = !1T 


(3-1) 


the skew symmetric matrix Q. being given by 


“3 -UJ 2 


Q - 


-U)r 


w 2 ' W 1 


(3-2) 


A trivial consequence of the equation expressing T in terms of T and the components 
of <u obtained by a simple transposition of its left and right members, is the 
relation 
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• T T T 

t = t a 


(3-3) 


which is used repeatedly in the development of the equations of motion. It is worthy 

T 

of note in passing that the skew symmetric matrix has the property that 


T -*■ 


B 


1 

o 

1 

8 

CO 

8 

to 

1 


" V B 1 ' 

COg 0 -03^ 


V B2 

. - ^2 0 _ 


- V B3 -1 




(3-4) 


—T —r — r — r 

where V' B = Vg^i + Vg 2 j + Vg^k AS an ar bi trar y vector expressed on the B (or B) 
vector basis. 


With flexible appendage i (1 < i < NA) one associates the following rotation 
matrices . 


T. = T(B -> i) q _ q is a prescribed rotation matrix of constants, peculiar to a specific 

vehicle configuration, defining the transformation from the B (or B) resolution 
to the i resolution when 6. = 0, that is, when the i-frame is in its "null'' orien- 
tation relative to the B -frame (or B -frame). 


1 = 
l 


[6 i ] (l) 



0 0 

cos 0. sin 0. 

-sin 0j cos 0j - 


(3-5) 


T. = J/T. = T(B-*i), the rotation matrix defining the transformation from the B (or B) 
resolution to the instantaneous i resolution. It is not difficult to establish the 
useful relations 







a. T ) 


(3-6) 


where 



ro 



0 

0 

1. 


0 

-1. 

0 


(3-7) 
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It is easy to see further that with the notation ub = 0,-h, the symbol ub denoting the 
angular velocity of the i- frame relative to the B- or B- frame (and expressed on the 
i vector basis) , 



V. = 
1 


x V. 
1 


(3-8) 


for any vector V. expressed on the i vector basis. 

+ Vi 

Pertinent to the i SDOF CMG are the following transformations, definitions, 
and consequential relations. 


T g ^ o = T(B-J-Gio), a prescribed rotation matrix of constants defining the transforma- 
tion from the B (or B) resolution to the Gi resolution when the gimbal angle 


6 Gi= °' 


The direction of K . (and hence of K~.) when 6„. = 0 and that of 
gi Gi Gi 

(Here, by direction is meant direction 


i„. determine the transformation T„. . 
Gi ^ Gio 

in the B- (or B) frames.) 



' 1 . 

0 

0 

*Gi [<S Gi ] (l) 

0 

cos 6 Gi 

sin 6 Gi 


. 0 

-sm « Gi 

cos 6 Gi . 


T G . = T g .Tq. q = T (B->Gi), the transformation from the B (or B) resolution to the 
Gi resolution. 

T . = = T(GVgi), the transformation from the Gi resolution to the gi resolu- 

tion. 

T . = T .Tq. = T(B+gi), the transformation from the B (or B) resolution to the gi 

o o 

resolution . 

With the definitions: 

> » -> ~ 

^Gi = ^Gi *Gi = an S u l ar velocity of Gi frame (gimbal frame) relative to the B and B 

frames expressed on the Gi vector basis (3-10) 


fi Gi = «Gi 


0 

0 

Lo 


0 0 

0 - 1 . 

1 . 0 


(3-11) 
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( 3 - 12 ) 


to • = <p . K . = to . K . = angular velocity of the gyro element (gi frame) 
& o 1 & 1 & g relative to the Gi frame 


ro -i. 


a 



l. 

o 


o 

0 


0 

0 

0 


there follows. 


• T 
T 

Gi 

T 

= T . 
Gi 

T 

fi Gi 




T 

Q Gi 

^Gi = 

“Gi X 

^Gi 

(for any vector 
vector basis . ) 

— > 

Vqj expressed on the Gi 

T 
Q . 
gl 

\i ” 

to . x 
gi 

V . 
gi 

(for any vector 
vector basis.) 

V . expressed on the gi 

4 T 
T . 
gi 

~ T 
= T . 
gi 

T 

fi . 

gi 




• T 
T 

gi 

T 

= T . 
Gi 

T 

S Gi 

~ T 
T . 
gl 

T T 

+ x . n . 
gi gi 

• 


Among other SDOF parameters presenting themselves in the derivation o 
moment equation and gimbal equation are 


,Gi 


I 


Gi 

"xx 

0 


0 


0 

I Gi 

yy 

o 


i 


o 

o 

Gi 
'zz J 


= inertia matrix of gimbal referred to y^j z^ 


r^ 1 = 


I ^ 1 

xx 

0 

0 


I 


0 

gi 

yy 

o 


o 

0 

zz 


inertia matrix of gyro element referred to 

x . y . z . (I^ 1 = I^ 1 assumed) 
gi J gi gi xx yy 


an important consequence of the assumption 1^ = 1^ 


being 


( 3 - 13 ) 

( 3 - 14 ) 

( 3 - 15 ) 

( 3 - 16 ) 

( 3 - 17 ) 
( 3 - 18 ) 
f the 

( 3 - 19 ) 

( 3 - 20 ) 
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(3-21) 


T . I gl T . = I gl 


gi 


gi 


from which it follows that 


T . I 8 ' 1 T • = T^. I gl T„ . 


T 

gi 


gi Gi 


Gi 


(3-22) 


Also appearing in the moment equation are the vectors w-,.' , Wq.’ , and w , the 
prime indicating the B (or B) resolution. They are defined by 


T t — T X T 

u,-,.' — T = T„ . T„ . . = T ~ . Up. = Tp . 

Gi Gi Gi Gio Gi Gx Gio Gi Gio 


“Gi’ = T Gio 


6 Gi 

0 

0 


5 Gi 

0 

0 


T x ~ T T 

w •' = T . w . = T„. T . u . = u . T 0 . 

gl gi gi Gi gi gl gi Gio 


-sin 6 g . 
cos 6^ 


(3-23) 


(3-24) 


(3-25) 


Transformations and subordinate relations pertinent to a 2DOF CMG follow. 

TgGB = T(B-GB), a prescribed rotation matrix of constants defining the transforma- 
tion from the B (or B) resolution of a vector to its gimbal base (GB) resolution. 

Tgbqg = T(GB^OG) = f^oG^(3)’ transformation from the GB resolution to the 
OG resolution. 


Tqgig = T(OG^IG) = t ^ jq 1 ( i) » the transformation from the OG resolution to the IG 
resolution . 


Tj G g = T(IG->g) = [ sPg. 5(2)’ transformation from the IG resolution to the g 

resolution . 

T BOG = T t B ^ 0G) = T GBOG T BGB = [lS OG^(3) T BGB = T ( B ' >OG )’ the transformation 
from the B (or B) resolution to the OG resolution. 
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T BIG T(B+IG) T(Sh-IG) T OGIG T GBOG T BGB [<S IG ] (1) [6 OG ] (3) t bgb ’ the 
transformation from the B (or B) resolution to the IG resolution. 

Tgg = T(B^-g) = TjQgTgjQ, the transformation from the B (or B) resolution to the 
g resolution. 

With the symbols Wqq , “jq> and w denoting the angular velocity of the OG- 

frame relative to the GB -frame (and hence relative to the B- and B-frames), the 
angular velocity of the IG-frame relative to the OG-frame, and the angular velocity 
of the g- frame relative to the IG-frame, respectively, one has 

^ . -> 

w OG = 6 OG K OG 

W IG = ^ IG X IG = ^IG X OG (3-26) 

w g = “g \ = w g j IG 

and the associated skew symmetric matrices 

0 - 1 . 

1. 0 

0 0 

0 0 
0 0 
0 1. 

0 0 1 . 

0 0 0 

-1. 0 0 

having the properties 

^OG V OG = t)G X V OG ^ f ° r arbitrar y v qg ex P ressed on the OG-vector basis), 
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ftj g V IG - i^ G x V (for arbitrary Vj G expressed on the IG-vector basis), 


V x V (for arbitrary V expressed on the g-vector basis). 

& & & £■> ^ 


From the definitions above one can easily deduce that 


• T T T 

T L = T o A 

1 GBOG X GBOG OG 


• T T T 

T = T o 

1 BOG 1 BOG ^OG 


■ t T T 

T = T o 

1 OGIG 1 OGIG “i G 


. rp rp rp rp rp 

T 1 = T 1 (T 0 x T 1 + 0 A } 

BIG BIG U OGIG “OG OGIG 


( 3 - 28 ) 


T ?Gg T IGg ^g 

*2g = < T BOG “Sg T OGIG + T S° G T ° GIG ^ G> T ‘ G « + T ® IG T ^ G « ’ 


these relations being of considerable utility in subsequent derivations. Regarding the 
derivations to follow, it should be remarked that the inertia matrices of outer gimbal, 

OG IG g 

inner gimbal, and gyro element, denoted, respectively, by I , I , and I s , are 

supposed diagonal with 1^ = 1^^ , 1^ = 1^, and 1^ = 1^ . The matrices , 1^, 
" & xx yy yy zz xx zz 

and I g are, incidentally, referred to x 0 G y OG z OG’ x IG y IG Z IG’ and X g y g Z g’ res P ec ~ 

tively . Abbreviations of certain combinations of these matrices appear in the moment 
equation. They are defined as follows. 


jOG+IG+g rpT jOG rp rpT / T _j_ T S? \ rp 

B BOG 1 BOG A BIG U 1 ; BIG 


lIG+g = jIG + jg 


jl G+g _ rp T jIG+g rp 

X B BIG BIG 


( 3 - 28 . 1 ) 


T& - rpT .g rp _ rpT rg rp 

B 1 Bg 1 A Bg BIG 1 1 BIG 
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The trace of 1^ is denoted by Tr(I^), that is, 


Tr(I ? ) = 1 K + + I ? , § = OG, IG, g 

xx yy zz ^ & 


(3-28.2) 


The first time derivatives of the vectors in equation ( 3- 26) , as measured in and 
expressed on the indicated vector bases, are obviously 


(dt) 


OG 


w OG w OG 6 OG K 0G 


(*) 


IG 


W IG W IG 6 IG \G 


(3-29) 


(a d t) 


g 


“g = W g = W g ] g = ° (u g = constant) 


The B (or B) resolutions of the vectors of equations' ( 3- 26) are given by 


rt * — t t 

OG BOG w OG 


w ig' T BIG w IG 


(3-30) 


w ' = tJ u 

g Bg g 


while the first time derivatives (as measured in the B- or B -frame) of the primed 
vectors of equations' (3-30) are given by 


(*) 

U) 


B 


B 


-* ! _ “*• , rpT -► 

w OG w OG x BOG OG 


T — , . T — 


W IG " W IG “ "BIG '"IG ' '"OG ~ '"IG 


1n T n X U). 


(3-31) 


(^) V = V = (a W + “IG^ x V 
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Attention is called here to a remark concerning derivatives that appears as the first 
sentence of the paragraph following equation (4-19). 

i. u 

Pertinent to the i in swivel engine is the transformation T^. = T(B->Ei), the 

transformation from the B (or B) resolution to the Ei resolution, whose time deriva- 
tive is shown in Reference 1 to be given by 


T„. = T„. - T„. q 

Ei El El Ei 


(3-32) 


where 


n 


Ei 



- o> 


(3) 

Ei 

( 2 ) 

J Ei 


- w 


(3) 

J Ei 

0 

( 1 ) 

Ei 


w 


( 2 ) 

J Ei 

( 1 ) 

Ei 

0 


(3-33) 


the to^P , j = 1, 2, 3, being the components of the vector which is here the 

angular velocity of the Ei-frame relative to the S-frame and is supposed expressed 
on the Ei vector basis, that is. 


Ei 


= T 


Ei (w B 


u E i') 


(3-34) 


the symbol cog.' denoting the angular velocity of the Ei-frame relative to B (or B) 
frame and expressed on the B (or B) vector basis. Among other relations satisfied 


-A. 

• 

by are the following. 




^Ei T Ei Q Ei' 




(3-35) 

^Ei = ^Ei T Ei 





where 

' 0 

'(3) 

“Ei 

'(2) I 

“Ei 


/ 'T 

fi Ei = - Q Ei = 

-n ,(3) 

“Ei 

0 

,;(i) 

“Ei 

, (3-36) 


/-v 

w W 

3 

I 

'CD 

“Ei 

0 
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and 



I 

O 

- (3) 
“Ei 

-hi* 1 


°Ei “ ~ Q Ei ~ 

-n (3) 
“Ei 

0 

hi* 

, (3-37) 


M 

“Ei 

-M (1) 

“Ei 

0 



the ^ , j = 1, 2, 3, being the components of and the . j = 1> 2, 3, those 

of the vector tJg. = T R . w R '. which is the angular velocity of the Ei- frame relative to 
the B (or B) frame expressed on the Ei- vector basis. An interesting consequence of 
the last two expressions for T R . is the similarity transformation, 



Q Ei T Ei 


(3-38) 


Complete specification of T R . and cannot be made until a "B to Ei Euler sequence" 

is prescribed, that sequence being dependent upon the point of application and direc- 
tion of the thrust delivered by the engine, the engine actuator arrangement, the sign 
convention for engine deflections, etc. 

Ei 

The symbol I appearing in the moment equation denotes the inertia matrix of 
swivel engine i referred to the axes x E iY E i z E j (defined in section 2). 

Transformations pertaining to the i n rotor are the following: 


T R . = T(B-*Ri), the rotation matrix of constants defining the transformation from the 
B (or B) resolution to the Ri resolution. Its elements are functions of the 
components of the unit vector A R . e [X p .,y n .,Vp.] which specifies the posi- 


tive direction (relative to the B-frame) of the spin axis (the x„. axis) of rotor 

Jtvl 


1 . 


A frame oriented as the B-frame could be carried into the Ri orientation by 

the 3, 2 sequence through th e angles 8 R . = tan > 0 < < 2 tt, and 

- 1, I /* 2 % 


Ri where ^ R . = tan 


-ip. 

determined by 




Ri 


'), -tr/2 < tjj R j ^ tt/2, in which case T R ^ is 


^Ri t _ 'f'p n'l I' 


f Ri ( 2) 1 Ri (3) * 


(3-39) 


= T(Rf>(Ri)*) = [v? R .](i), the transformation from the Ri to the (Ri)' resolution. 
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T„. =3-0- T„. = T(B->-(Ri)') , the transformation from the B resolution to the (Ri)' 
Kl Kl Kl 

resolution . 


From the above it is a trivial matter to show that 


Ri i Ri u Ri 


^ rp ^ rp _rp ^ 

T Ri = T Ri Ri + n Ri ) ’ 


The skew symmetric matrix being given by 


Ri Ri 


0 0 0 

0 0 - 1 . 

0 1.0 


w Ri ^Ri 


with the obvious property that 


( 3 - 40 ) 


( 3 - 41 ) 


~T 

o V 
Ri (Ri)’ 



V 

(Ri) T 


( 3 - 42 ) 


where the vector V^.^ t is an arbitrary vector expressed on the (Ri)' vector basis 
and denotes the angular velocity of the (Ri)'- frame relative to the B- frame 
expressed on the (Ri)' vector basis. By definition of one can write 


OJ 


Ri 






( 3 - 43 ) 


the scaler v 3 ^. = admitting of both positive and negative values to correspond to 

positive or negative rotations (in accordance with the right-hand rule) about the rotor 
spin axis. Obviously, 


U) 


Ri 







( 3 - 44 ) 


and furthermore, if denotes the B (or B) resolution of the angular velocity of 

the (Ri)' frame relative to the B (or B) frame then 


(0 


Ri 


T 



Ri 


( 3 - 45 ) 
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and it follows that 


Ri 


■(») 


B 


> T —> ~ r T^ -> m -4 

“Ri' = T Ri n Ri u Ri + T Ri ^Ri = T Ri w Ri = T Ri w Ri 


(3-46) 


Ri 

The inertia matrix, I , of rotor i referred to x' -y' .z' . is herein assumed 
Ri Ri RI RI rll 

diagonal with I = I zz - Under these assumptions it follows that 


7 T . i r W = i Ri 

-'Ri J Ri 1 


(3-47) 


and hence that 


rpT ,Ri px rpT T Ri rt ty _ rpT .Ri 

Ri Pi 1 Pi“ / T3i 1 -'Oi 1 T5 “ 1 o ,■ 1 1 D4 > 


Ri Ri Ri 


'Ri Ri 


Ri 


Ri 


(3-48) 


a fact used in arriving at equation (4-67) in the following section. 


Passage from the S resolution to the E resolution is effected via the transforma- 


tion 


Tgg T(S-*E) [ otp + uig (t t Q )] ^ , (3-49) 

the symbol otp denoting the right ascension of the prime meridian at time t Q and w e 
the magnitude of the Earth's spin vector. 

Not to be overlooked, and certainly not the least important of the several trans- 
formations discussed herein, is the relation between the unit vector triad u D u u r 

R A o 

associated with a point of the earth's exterior gravitational field and the triad f c ? c K , 

o o o 

that being 


1 

C i 

w 

1 


"V 




U A 

[-6](2) 

J s 

-X 


1 

Ul 

l 


where a and 6 

point. If R = 
s 


are the right ascension and declination, respectively, of the field 
T 

[Xg, Yg, Zg] - Xglg + YgJg + ZgKg denotes the position vector of 
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the field point referred to the S-frame then the angles a and 6 are given by 


a = tan _1 (Y s /X s ) , 

6 = sin 1 (Zg/Rg) 

If A is the longitude (positive 
also given by 


0 <_ a < 2tt , 

-5 < 4 <1 - R s= |S s l 

east of the prime meridian) 


= A g 2 + Yg 2 + Zg 2 (3-51) 

of the field point then a is 


a = ap + oo e (t-t Q ) + X . (3-52) 

The unit vectors Ug, u-^, and u^ can also be expressed in terms of Ig, Jg , and Kg 
as follows. 

U R = V R S 
- K g x R g 

U, = , R q 4 ±R~ K , (3-53) 

A |Kg x Rg| b b b 

u 6 = u R x u A . 


Relative to the rotating E-frame (the earth fixed frame X-Y.-.Z-) the field point has 
-*■ T E E E 

position R £ = [Xg, Yg, Z £ ] = T(S->E) Rg so that A and <S can also be computed in 

accordance with 


A = tan 1 (Yg/Xg) , 0 < A < 2tt 


(3-54) 


6 = sin 1 (Zg/R) 


IT . . IT 

2 - 6 < 2 


R - Rg| - Rg 


and in terms of Ig, Jg, and Kg the unit vectors Ug , u-^, and u fi are given by 


_ ! r ' 


*E 

U A 

t-S]( 2 ) [^] (3) 

J E 

-* 


-* 

U o 


K, 

L 6 


E_ 


(3-55) 
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In Reference 6, the acceleration due to the Earth’s gravity at the point with spherical 

coordinates (R,A,<5) referred to the Earth fixed frame X^YgZg is resolved into its 

radial, longitudinal, and latitudinal components, these being the u^ component, the 
— ► — ► 

u, component, and the u A component, respectively. 

A 0 

SECTION IV. THE EQUATIONS OF MOTION 
The coordinates defining system configuration include the following: The com- 

S ~ ~ ~ -T> 

ponents of the vector R = [X , Y , Z ] , the position, referred to the S-frame, of 

s s s s 

the origin of the B -frame (the structural axes xy z) ; the Euler angles and 'P ^ 

specifying the orientation of the B -frame (and also that of the B -frame) relative to 

the S-frame; the angles 6., i = 1, ... NA, the symbol 6, denoting the angle through 

which the l flexible appendage* is permitted to rotate as a ” whole" relative to the 

rigid central carrier; the generalized bending displacement coordinates rw, j = 1, ... 

th 

N., associated with the i flexible appendage, i = 1, ... NA; the displacement 5 ., 

i = 1, ... NP, of the "point" mass m . from its equilibrium position; the angular 

fh ^ , 

deflections g . and g . of the i swivel engine, i = 1, ... NSE; the gimbal angle 6^. 
tvi y^ pA 

of the i SDOF CMG , i = 1, ... NSDOF ; the outer gimbal angle 5 q G . and inner 
gimbal angle 6 of the i* R two DOF CMG, i, ... N2DOF. 

The principal of virtual work finds itself of considerable utility in developing 
the equations of motion. That principle**, as applied to dynamic systems, states that 
in an arbitrary virtual displacement of the system (compatible with the constraints) 
the virtual work done by the inertial forces plus that done by the external forces 
equals the change in strain energy. On invoking the principle of virtual work and 
appealing to the independence of the generalized coordinates defined above, it follows 
that 



3 (MOON) 
R S 

(MOON) | 3 


+ 


R 


R (SUN) 


(SUN) | 3 
S 1 


dm + Qq 


3D 

3q 


3 U 

3 q 


(4-1) 


where q may be any of the aforementioned coordinates. 


* Later in this section, attention will be directed to a system wherein each flexible 
appendage is allowed two rotational degrees of freedom. 

** See Reference 16, pp . 114-115. 
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In the volume integral in the left member of equation (4-1), R denotes the 

s *• 

position referred to the S- frame of a generic point of the vehicle, while 1$ denotes 

s 

the second time derivative, as measured by an observer in the S -frame, of the vector 

R . (It should be remarked here that the number of dots above a vector symbol 
s 

indicates the order of the time derivative of that vector as measured by an observer 

— ► 

in the reference frame on which the vector is expressed.) The superscript T on R n 

_> np -> ® 

in equation (4-1) indicates that R is the transpose of the 3x1 column vector R . 

s s 

It is to be understood that when used as a subscript on an integral sign as in equa- 
tion (4-1) above, the letter m simply indicates that the integration extends over the 
volume occupied by the entire vehicle system, while if it appears as one factor of a 
product, it denotes the numerical value of the mass of the entire system. 

The generalized force is such that the product Q^6q is the virtual work 
done by those forces not derivable from either the potential* function U ( as herein 
defined ) or the dissipation function D when only the coordinate q undergoes the 
virtual displacement 5q, the forces alluded to being ( in this paper ) the gravitational 
forces of Earth, Moon, and Sun; direct solar radiation; aerodynamic forces; the 
torques applied to rotate the flexible appendages (relative to the rigid central carrier) ; 
and the torques applied to rotate the CMG gimbals. Accordingly, Q must be given 
by 


f - ( 9 *s 

Qq ~ (m) (d Gravity, + ~Jq~ (d ^ero.S'* 


f (A V x + // 

(A) 3q (d F solar, S'* 


(4-1A) 


where 


F GRAVITY, S 


^(EARTH) 
GRAVITY, S 


£(MOON) 
GRAVITY, S 


S(SUN) 
r GRAVITY, S 


= GRAVITATIONAL FORCE EXERTED BY EARTH, MOON, 
AND SUN ON THE ENTIRE SYSTEM 


*The symbol U in equation (4-1) is better termed the "strain" energy function. Use 
of the term "potential function" should not mislead the reader to believe that "every" 
force for which there exists a potential function "should" be derivable from U. 
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= AERODYNAMIC FORCE ON THE ENTIRE SYSTEM 


F AERO, S 

FgoLAR S = RESULTANT 0F DIRECT SOLAR RADIATION PRESSURE ON 
SURFACE AREA 


the subscript S again being indicative of the S-resolution, and where, necessarily , 


= 0 if< 


(4- IB) 


q # 0. , i = 1, ... NA 
q * S G ., i = 1, ... NSDOF 


q * <5 0Gi » i = 1, ... N2DOF 
v q * S IGi , i = 1, ... N2DOF 


the symbol Ji denoting the torque applied to impart a change in the coordinate q 

q 

when q is any of the angular coordinates 0 ., 6q., Sqq., or <5 jqj (for each i over the 
respective ranges) . The significance of the subscript m on the first integral sign in 
the right member of (4-1A) has already been discussed. It should be self evident 
that the subscript A on the second and third integral signs in (4-1A) simply indicates 
that the integration extends over the surface area of the entire vehicle (actually that 
portion of the surface experiencing impingement of air molecules in the case of the 
second integral and that part of the surface exposed to direct solar radiation in the 
case of the third integral) . 


The most simple of all approximations to Fq^aVTjY g 
will here be considered satisfactory. They are 


and F< SUN) 

r GRAVITY, S 


S(MOON) ~ 

GRAVITY, S my m 


-(MOON) _ j^CM. . i — (MOON) _ £CM,3 
^ s s } ' S S I 


(4- 1C) 


S(SUN) 
GRAVITY, S 


-(SUN) 
my (R 

M s s 


g CM |g (SON) _ gciy 

S ' 1 S S 1 


(4- ID) 


where 


y m = product of universal gravitational constant and mass of Moon 

y = product of universal gravitational constant and mass of Sun 
s 
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I* (MOON) = position vector of moon referred, to the S- frame 

R^SUN) _ p 0s ^i on vector of sun referred to the S-frame 

■pCM _ S . rpT ~ 

R s " R s + T r CM 


?cm = position of vehicle CM referred to the B- frame, 


the components of r( moon ) an( j ^(SUN) be ing presumed known tabular functions of 

S S 

time, that is, available from the ephemerides of moon and sun. The far more detailed 

-'Y'F ARTFiH 

approximation to the B -resolution of the vector ^ GRAVITY g> herein denoted by 
Fgg , is provided by Appendix A . 

The vector sum 


M 


R 


(MOON) 


Pc R 


(SUN) 


iS(MOON) 3 

| R S 


;( SUN) | 3 


(4- IE) 


appearing within braces in the left member of equation (4-1) is merely an approxima- 
tion to the acceleration of the S-frame relative to the Newtonian frame (the N-frame) 
and is realized by treating Earth, Moon, and Sun as point masses and ignoring all 
forces on the earth other than the gravitational forces of moon and sun. See Figure 1. 

Subtracting m times the first term of the expression (4- IE) from both members 
of (4- 1C) followed by a pre- multiplication by the rotation matrix T = T (s-*B) gives 
the vector 


, g (MOON) gCM ^(MOON) 

it (MOON) _ m T I_S ~ S __S 

gB U M 1 1 .-(MOON) S CM I 3 | ;t (MOON) i3 

l |R S - R g | |R S I 

appearing in the translational equation below. A similar combination of the second 
term in expression (4- IE) with the members of (4- ID) gives, after a pre-multiplication 
by T , the vector 

/ k(SUN) gCM s(SUN) 
it(SUN) _ _ ^ ) S " K S K S 

e B “ s (|r< suN) - i?™| s ‘ jg<'s<W 




(4-1G) 




(4- IF) 
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also appearing in the translational equation. 


One will notice that regarding the S -frame as Newtonian (as the author has 
seen some people do) would result in the omission of the expression (4- IE). Although 
the S-frame does not rotate relative to the N-frame, it still cannot qualify as Newton- 
ian by virtue of its acceleration relative to the N-frame. 

The author of this paper has no comment regarding the significance of the 
error introduced by the complete neglect of electromagnetic forces; reflection of solar 
radiation by the earth and its atmosphere; direct thermal radiation from the earth; 
micrometeorite impacts, the attraction of celestial bodies other than the Earth, Moon, 
and Sun; etc. 


The thrust is embedded in the first volume integral and the generalized force 

Q in equation (4-1), the dominant part (often called the "momentum" component) 

4 

belonging to the volume integral and the other part (referred to as the "pressure" 
component) belonging to the surface integrals which comprise a part of Q . 


The major contribution of the i 
is Hie strain energy U^. given by 


,th 


flexible appendage to the potential function U 


N. 


U 


Ai 


.2 .2 

n/il 1 1 

Mr « k n K 


K=1 


i = 1, . . . NA 


(4-2) 


small deformation theory having been assumed and thermal effects ignored. The 
symbol oA denotes the frequency of undamped free vibration in the K th natural mode 
whose shape function, the 3x1 column , is a function of position r. = [x. , y., 
z,] 1 , referred to the i-frame (defined in section 2), of points of flexible appendage i 
in its undeformed state. The generalized mass, Mr, associated with the K l natural 
bending mode of flexible appendage i is defined by 



K = 1, ... N., i = 1, ... NA 


(4-3) 


A remark regarding the symbol m., similar to that regarding m, is in order here, 
that being that when m^ appears as a subscript on an integral sign as in equation 
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(4-3), it indicates that the integration* extends over the volume occupied by flexible 
appendage i, while elsewhere nn denotes the numerical value of the mass of flexible 
appendage i. 

An approximation to D^, here defined as one-half the rate at which energy is 
dissipated through structural damping in flexible appendage i, is 


N. 


D 


Al 


% 2 


£ 

K=1 


2 M 


K "K 


i = 1, . . . NA 


(4-4) 


• th 

the symbol ^ denoting the damping ratio associated with the K L natural bending 

mode of flexible appendage i. The approximation (4-4) is based upon the assumption 
of the existence of a viscous damping coefficient distribution (per unit volume) which 
varies directly as density (Reference 1) . Although it is generally agreed that the 
viscous approximation to structural damping is inadequate, the damping force deriv- 
able from equation (4-4) will herein be deemed satisfactory. On summing the expres- 
sion for D . over i (i = 1, ... NA), one has an approximation to the contribution of 
structural damping in all of the flexible appendages to the dissipation function D . 

An elastic restoring force -K Q . 0. and viscous damping force -C Q . 0., both of 
which oppose the rotational motion (relative to the central carrier) of flexible append- 
age i, give rise to the contributions 

U 0 . = (1/2)K 0 . 6. 2 , i = 1, ... NA , (4-5) 


and 


D ei=< 1/2 > c ei h 2 


i = 1, . . . NA 


(4-6) 


to U and D, respectively. 


The restraining spring and viscous damper forces, ~K jCpj and 


C pi^pi’ 


which 


impede the displacement (£ .) of the point mass m . account for the following terms 

pi pi 

in the expressions below for U and D . 


*Context will make clear whether the integration extends over the undeformed or 
deformed appendage. 
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V =a/2, V Si • I = 1 -"- NP - 


(4-7) 


V = U/2) C pi «pi ’ i=1 --- NP 


(4-8) 


To account for the energy loss through viscous friction in the gimbal bearings 
of the CMG's, one should include in the expression for D terms pertinent to the SDOF 
and 2DOF CMG’s, these being, respectively, 


°6Gi = (1/2) C Qi <5^. , i = 1, ... NSDOF , 


(4-9) 


and 


D 6 OGi C 1 / 2 ) C QGi 6 OGi 


i = 1, ... N2DOF 


(4-10) 


D 6IGi =(1/2)C IGi hw ’ i = 1. ... N2DOF . 


(4-11) 


The terms in U and D pertinent to swivel engine i are, respectively, 


U SEi 


2 [K PEi (3 pi ” ^PCp + K YEi (3 Yi “ 3 YCi ) ] 


(4-12) 


and 


D SEi 


5 [C PEi »Pi + C YEi '40 ■ i=1 > ••• NSE 


(4-13) 


In equation (4-12), Kp^ and Ky^ are "effective" torsional spring constants having 
the dimension NEWTON *METERS/RADIAN . The Cpgj and Cy^,. equation (4-13) are 
"effective" viscous damping coefficients (in pitch and yaw, respectively) having the 
dimension NEWTON *METERS/ (RAD /SEC ) , that of a viscous torsional damping coefficient. 
The angles Bp^ and are engine deflection commands (in pitch and yaw, respec- 

tively). From equation (4-12) is derivable the moment about the swivel point imparted 
by the actuators, and from equation (4-13) the damping moment retarding engine 
deflection. From equation (4-2) and equations (4-4) through (4-13), there follows 
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(4-14 


+ 2 ^ / ^PEi^ ^Pi ^pci^ + 3v/->i) 1 


YEi Yi YCi 


i=l 


D 


NA 

NA 

NP 

NSDOF 

N2DOF 

E d m 

i=l 

+ E d 6. 

i=i 1 

+ E D Pi + 

i=l 

e 

i=l 

D. + YI + > 

6 Gi . t 6 OGi 6 IGi 

1=1 


+ 


1 

2 


NSE 



E D SEi 



i=l 



XT A N. 

NA l 

9 NA 

NP 

E E 2 < 

i i • i , 1 V r 
C K W K n K 2 2^ C e. 

°ME C 

i=l K=1 

i=l 1 

i=l 


+ 


NSDOF N2DOF 

2 E C Gi ^Gi + 2 E (C OGi ^OGi + C IGi ^IGp 
i=l i=l 


NSE 


+ 2 (C PEi ^Pi + C YEi P Yi ) 


i=l 


(4-15 


In developing the equation of motion corresponding to a system coordinate q in 
accordance with equation (4-1), it is convenient, if not necessary, to first write the 
volume integral /* (3R c /3q) T R„ dm as the sum of several integrals. 



f 


-►T 

9R S 

8q 



dm 




NSE 


NP NA 




(/ 

♦ f 

* Z 

/ + 

Z / +Z 

/ 



\<”o> 

(m f ) 

i=l 

(m Ei ) 

i=l (m pi ) i=l 

(m.) 




NR 


NSDOF 

NSDOF 


N2DOF 


+ Z 

/ ♦ 

Z 

/ * z 

/ 

♦ z 

f 


i=l 

(m Ri ) 

i=l 

<m Gi ) i=l 

(m' .) 
gi 

i-1 

(n, OGi> 


N2DOF 

/ 

N2DOF 

* z 

+ z 

i=l 

( m lGi> 

i=l 



(m .) 

gi 



3q 



dm 


(4-16) 


each pertinent to a particular part of the vehicle and sharing with that part an 
intimate connection with one or more system coordinates. Equation (4-16), wherein 
the notation is self-explanatory, is complemented by the following definitions: 


m Q = mass of the rigid central carrier 

m f = mass of liquid propellant (in the tanks and feedlines) plus the products 
of combustion in all the engines forward of their respective nozzle exit 
planes (the word "forward" as used here is that of the thrust delivered 
by an engine) . 

t h 

nig. = mass of the 1 swivel engine. 

+ V\ 

nip. = mass of the 1 "point" mass. 

x l— 

m. = mass of the 1 flexible appendage. 

m Ri ~ mass of the i th rotor (the i th rotor may or may not be a reaction wheel). 

m^. = mass of gimbal of i th SDOF CMG . 

m gi = mass °f gyro element of i^ SDOF CMG . 

m OGi = mass °f outer gimbal of i th 2 DOF CMG. 

m IGi = mass °f inner gimbal of i in 2 DOF CMG. 

m gi = mass °f gyro element of i^ 2 DOF CMG . 
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The manipulations leading to the system equations of motion begin with the 
expressions for R and its requisite derivatives in terms of the generalized coordi- 

o 

nates and their time derivatives, in general, the position of an arbitrary point of 
the vehicle relative to the S- frame (defined in section 2) is given by 

H s = § s + tT ? ’ (4 ~ 17 > 



the vector r being the position of the point referred to the B- frame (section 2) and 
T i T (S+B) = T (S+B) the rotation matrix defined in section 3. Direct differentia- 
tion then gives 

R = R + T T (r + n T r) (4-18) 

s s 

.. * o 

R = R + T 1 [r + 2 r + (fi 1 + P 1 ) rj , (4-19) 

s s 

use having been made of equation (3-3) and the skew symmetric matrix defined by 
(3-2) with the property (3-4). 

To be recalled at this point is a remark made earlier to the effect that the 
number of dots above a vector symbol is indicative of the order of the time derivative 
of that vector as measured by an observer in the reference frame on whose associated 
vector basis the vector is resolved. Thus, 



In view of this convention and the property (3-4) of q , it should be evident that the 
expression within braces in the right member of equation (4-19) is completely equiva- 
lent to the more familiar "textbook like" expression 
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(i)rd\ 


r + 2 Wg x 


(d d t)g ? + [(dt)g 3 b] 


x r + o)g x (oig x r) 


(4-20) 


all vectors in equation (4-20) having the B -resolution. Thus, with the "derivative 
convention" of this report, one can dispense with the cumbersome notation of equation 
(4-20) and pass from equation (4-17) to equation (4-18) and from equation (4-18) to 
equation (4-19) as if in contempt of (yet in complete accord with) the differential 
operator relation 

(ar) s = (ar) g + “b x • 


Clearly, the only system coordinates (and derivatives thereof) in evidence in 

equations (4-17), (4-18), and (4-19) are the components of R and the Euler angles 

T T S 

(the Euler angles entering through T and ft which are completely determined once 
an Euler sequence has been prescribed). Consequently, it is necessary to comple- 
ment those equations by expressions for r, r, and r in terms of the other system 
coordinates and their time derivatives. Each expression alluded to pertains to a 
particular subdivision of the vehicle system, the manner of "subdividing" the system 
being dictated by the subscripts in equation (4-16). In what follows, context should 
make clear the intended meaning of the "mass symbols" m , m f , m„., ... m .. For 

O I ill g“l 

example, the symbol m Q in the paragraph immediately following is an abbreviation for 
the "rigid central carrier." 

At points of m Q , one has simply 

— ► — ► 

r = r 

~ ^ -*■ 

r = r = 0 , (4-21) 

• «• 

the vanishing of r and r being a consequence of the assumed rigidity of m Q and the 

definition of the B -frame as one at rest relative to m . 

o 

— ► 

At points of the fluid mass m^., the vector r is a function of both position and 
time, that is, 
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x (x, y, z, t) 


.x (x, y, z, t) 


so that 


(Vx) • r + 9 x/9 t 
(Vy) • r + 9 y / 9 t 


(Vz) • r + 9z/9t 


where 


V = i 9 / 9 x +j 9 / 9 y + k 9/9z 


Via the relations 


9f_ 9 , f 3p 

p w (pf) ‘ f It 


- V • (pr) (the continuity equation) 


V • (fp r) = f v • (pr) + (Vf) • (pr) 


where f is a scalar function of position (r) and time t, and p is density (also a func- 
tion of position and time) , one can go further to show that 


V • (x p r) 


(pr) + v • (y p r) 

• -* 

. V • (z p r) 


an important result (intimately related to thrust) finding its application in the deriva- 
tion of the translational and rotational equations. 


At points of m^,. 
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r r Ei " £ Ei A Ei + T Ei r Ei 



(4-25) 


i _ j- • T -► / 

r “ "*Ei A Ei + T Ei r Ei ’ ( 

®) H ?Ei 5 S 

9 

(4-26) 


11 •• 

~ -*• • ■ T 

r “ '^Ei A Ei + T Ei r Ei 



(4-27) 

-+ 






r Ei being the position (referred to the B -frame) of the swivel point, the positive 
scalar the distance between the swivel point and the engine CM; the unit vector 
Agj (having the B- resolution) directed as the thrust the position of a generic 

point of the engine referred to x £ . Zg. (defined in section 2); and T £ . the rota- 
tion matrix whose definition appears among those in section 3. 


Pertinent to the point mass m one has 

pi 


r = r . + £ .A . 

pi ‘’pi pi 


(4-28) 


i . A . 

Pi Pi 


(4-29) 




(4-30) 


r . being the equilibrium position of m . and £ , its displacement in the direction of 
pi pl pi 

the unit vector Ap^ if £p. > 0 and in the direction of -Ap. it < 0. The unit 

vector A . has a constant direction in the B -frame. (If the mass m . is a part of the 
pi pi 

mechanical analog of a consumable liquid, liquid propellant for example, then as the 

liquid is depleted, the vector r . will vary which is to say that its time derivatives 

Pi 4 ^ * 4 . 

will not vanish. However, even in such a case the approximations F % 0, r n . % 0 

will be made herein.) 


Pi 


Pi 


For dm € m., that is, at points of flexible appendage i, 

? = [?j + A (?. ; t)] (4-31) 

r = T? 1 [r. + A (?j, t)] + T? A(?j, t) (4-32) 
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(4-33) 


= T T 
1 


[r i + 


A(r., t)] + 


2 Tj 1 


(r., t) + 


(r., t) 


In equation (4-31) r. denotes the position of the point of attachment of flexible append- 
age i; the matrix T\ is as defined in section 3; r^ denotes position relative to the 
i-frame (section 2) when the appendage is in its undeformed state; and A (r^, t) is 
the displacement, due to deformation of the appendage, experienced by the point 
whose position (referred to the i-frame) before deformation was r.. The generalized 
bending displacement coordinates nj, j = 1, ... N^, are to be so determined that 
A (r.,t) is given by 


N. 


(?i. t) = L nj ?.. (i) 


3=1 


* (i) ^ (i > 


(4-34) 


the 3 x N. modal matrix $ ^ 
elements are functions of r^, 


having for its j 


•th 


and the N. x 1 column n having 


column the 

(i) 


3x1 column v 


(i) 


i Uh 

n j as its j 


whose 
element . 


Denoting the position (referred to the B -frame) of the CM of rotor i by r^. and 
position referred to the (Ri)*-frame (section 2) by r^, the position (relative to the 
B -frame) of a differential element of mass belonging to m^ must be given by 





(4-35) 


from which 






(4-36) 


it being presumed that r R; . is a constant vector on the B basis and that the rotor is 
rigid so that (d/dt)^. R .^ t r^d e 0. The "rotors" of this paper may or may not 
include reaction wheels but do not include the rotating elements of a CMG . 


The position of a differential mass dm t m' ., referred to the B -frame, is 

(jrl 


r = r Gi + T Gi r Gi 


(4-37) 
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where r^,. is the position of the CM of mj-,., r^,. denotes position relative to the 
Gi-frame (section 2), and T_. the rotation matrix whose definition (section 3) is 
found among those pertinent to the i SDOF CMG. With r Q j regarded as a constant 
vector (in the B -frame) and m£,. regarded rigid, there follows 



Similarly, for dm e 


m 


gi 



+ 







(4-39) 

(4-40) 


F; being the position of the CM of m' . (also assumed rigid) ; r . the position referred 
to the gi-frame (section 2); and T . as defined in section 3. (The assumption that 


m' . and m' 


have a common CM implies, of course, that r . = r„.. 

gi Gi 


Gi — "gi 

further that the gimbal axis passes through the CM.) 


It is assumed 


In the definitions and consequent relations pertaining to a two DOF CMG 
(Sections 2 and 3) the additional subscript i (suggestive of the i^ two DOF CMG) 
was omitted on all symbols incident thereto not only for "convenience" but deemed 
unnecessary by the author in view of the self-evident fact that the orientation of the 
associated coordinate systems and the structure of the relevant transformation matrices 
change from CMG to CMG . In keeping with the aforementioned definitions and rela- 
tions, the subscript i will be suppressed on all symbols in this paragraph, it being 
understood that they apply to a two DOF CMG. Thus, under the assumption that the 

outer gimbal, inner gimbal, and gyro element of the two DOF CMG in question are 

. . ^ ^ ^ ~ 
rigid with a common CM located at r„,_, = r T ~ = r relative to the B -frame, there 

Uu Hj g 

follows: 


dm € m 0G -* 


— ► 


„T 

r = 

r OG + 

T BOG r OG 

— > 

• T 

-► 

r = 

T 

BOG 

r OG 


■* T 

— ► 

r = 

T 

BOG 

r OG 


(4-41) 


33 


r r IG + T BIG r IG 


dm e m 


IG 


r — T v 

BIG r IG 


?= T T ? 
BIG r IG 


r = r + T„ r 
g Bg g 


dm e m -* 

g 


r = T F r 
Bg x g 


?=Tj ? 

Bg g 


(4-42) 


(4-43) 


where r denotes position relative to the 5 -frame, £ = OG , IG , g, and Tg^. = (TB-*0 
denotes the rotation matrix defining the transformation from the B (or B) resolution 
to the 5 -resolution , £ = OG, IG , g. 

Turning to the development of the equations of motion, attention will first be 

focused on the vector equivalent of the three scalar equations corresponding to the 

components X , Y , and Z of R . By mere inspection of equations (4-1), (4-17), 
s s s s 

(4-19), and the expressions for U, D, and Q (q = X , Y , Z ), it should be obvious 

Q S S S 

that the vector translational equation, before further simplification, reads as follows: 


f (R g + T T [r + 2 n T 
(m) 


m 


B (MOON) 
y m s 
g (MOON) ,3 


r + (0 ) r] } dm F GRAVITY S 

m p R (SUN) _ 

K s s ~ 

(SUN) | 3 + F AERO, S + F SOLAR, S 


(4-44) 


a result that one might have expected without resorting to equation (4-1). An imme- 
diate simplification of equation (4-44) is 


m T R = -m (fl^ + ) r 


CM 


' /< 


r + 2 q r) dm + F gB + F gB 


(SUN) 


(m) 


, £ (MOON) d? - 

gB F AERO F SOLAR 


(4-45) 
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All vectors in the right member of equation (4-45) have the B (and also the B) reso- 
lution. The vector F^ B is defined in appendix A, while and F^ B ^ UN ) 

are given by equations (4- IF) and (4-1G), respectively. Integral expressions for 
F AERO and F SOLAR are £* ven in appendices B and C, respectively. The vector 
? cm , the position of the CM of the entire system relative to the B -frame, is defined 
in appendix D. Via the relations (4-21) through (4-43), one can argue that the 
integrand of the integral in (4-45) vanishes except at points of m^, m^ (i = 1, ... 
NSE), m. (i = 1, ... NA), and at the point masses m^ (i = 1, ... NP). 

Integrating both members of equation (4-24) over V^, the volume occupied by 
m^, and applying the divergence theorem gives 



r dV 




^ (P r) dV 



r (p r • n f ) dA , 


(A f ) 


(4-46) 


where A^ is the area of the surface bounding and n^ is a unit vector normal to 
A^ (directed inward) at a generic point of A^. The integrand of the surface integral 
in equation (4-46) vanishes at all points of A^, except on A £ (the sum of the exit 
areas of all engines, both swiveled and fixed, associated with m^) . Thus, 



(4-47) 


the symbol n denoting the unit normal vector (directed inward) at a generic point of 
A„. The surface integral in equation (4-47) will be recognized as the "momentum" 
component of the thrust associated with m^. The thrust, F T , is defined by 



/ 


r(pr*n)dA + 


/ 


(P - P Q ) n dA 


(4-48) 


the symbols P and P q in equation (4-48) denoting, respectively, local static pressure 
and free stream static pressure. On dropping the rightmost surface integral (the 
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"pressure" component of F T ) in equation (4-48) as negligible and discarding also the 
first two terms in the right member of equation (4-47), the first being the Coriolis 
force and the second representing the unsteadiness of the fluid flow, one has the 
approximation 


/ 

(m f ) 


(r + 2 Q r) dm % F 


T ‘ 


(4-49) 


Manipulations (omitted here) will show that 


/ 

(m Ei ) 


(r + 2 Q r) dm = 


m Ei £ Ei ^ A Ei + 2 Q A eP ’ 1 


= 1, ... NSE, (4-50) 


/ 

(m Pi) 


(r + 2 r) dm = - m pi ( £ p . A p . + 2 £ pi A pi ) , i = 1, ... NP , 


(4-51) 


/<? 


(m.) 


+ 2 Q T r) dm = - rm T? {it x Y (l) n (l) ] 0. 


+ i. x [i^ x (&.^ + V ^ r/^)l 


+ 1 x [* (i) ; (i) ] 0,} 


- 2 m. q T tT {i! x [it/ 0 -* + e. 


+ - mj tT , i = 1, ... NA. 


(4-52) 


In arriving at equation (4-52), use was made of equations (3-6), (3-7), and (3-8). 
The symbol £. v J denotes the position, relative to x.y.z^ (the i- frame defined in 
section 2), of the CM of appendage i when = 0, that is, when the appendage i 

is in its un deformed state. The 3 x N. matrix 4* ^ has for its column the 3x1 
column matrix 
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= (l/m.) f dm 

/ *L \ 


> j~l> • • • N 


(4-52.1) 


C&j) 


which is to say, obviously, that 


¥ ^ = (l/nn) f dm 

1 (m.) 


(4-52.2) 


In view of equations (4-49) through (4-52), the translational equation (4-45) 
can be made to assume the form 

» t t. - - <« t * « t2 > ? cm ♦ ? gB ♦ > + ?<r oN> ♦ ♦ f aero 


NSE 

+ m £ . fcg. (A Ei + 2 0 A £j) 


SOLAR ' “‘Ei ^Ei v “Ei 

i=l 


NP 


E>i ( ^i r Pi + 2 ? 


pi a p1 ) 


i=l 


NA 


’ S m i ^ {i i X [ ^i (0) + ' i<(l) n(l)] ®i + \ x ^ (l)] ®i 


i=l 


+ i^ x [i. x + T n^)] 0.^} 


NA 


- 2 q T ^ m. t. t {i. x u/ 0) + v*'® n (l ^l e. + n^} 

1 1 1 X 1 


i=l 


NA 


- E m i l 


T *(i) ^(i) 


(4-53) 


i=l 


Note the absence of time derivatives of ? cm in equation (4-53). Therein lies the 
reason for choosing the components of R as three of the system coordinates instead 

S -* PM 

of the components of the position vector, R , of the CM of the entire system 

o 
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referred to the S-frame. Should one want the differential equation satisfied by 

-*■ fM Xj -* fM 'T' zt 

R , he has only to make the substitution R = R - T r in equation (4-53). 
s s s cm 

Definite expressions for and its time derivatives in terms of g p ^, and 
their time derivatives cannot be written until one specifies the actuator arrangement, 

— y 

the sign convention for g p ^ and gyj and the direction of the thrust when 3 p ^ and 
gy. are zero. If the use of swiveled engines proves to be impractical, as is likely to 
be the case, then, clearly, the terms reflecting their effect can be deleted without 
destroying the validity of the equations in which they appear. 

The effect of the rotors will not become apparent until the vector rotational 
equation has been written, this problem being addressed in the next and several 
succeeding paragraphs. 

The "customary" starting point in the approach to the moment equation is the 
generalized principle* of angular momentum (sometimes also called the generalized 
angular momentum equation) or special case thereof. Writing the moment equation via 
equation (4-1), as will be done in this report, requires no prior knowledge of the 
aforementioned principle; however, the use of equation (4-1) will require a few more 
manipulations, the extra preliminary manipulations serving in fact to establish a 
relation from which it is possible to deduce the principle of angular momentum (see 
Reference 8, pp . 20-21 of Chapter 2, paying particular attention to equations 14, 17, 
and 18) . 


Supposing the moment reference point** (MRP) to be the CM of the entire 
system, one should write (with r denoting position relative to the B-frame defined in 
section 2) 


-> “^CM T 

R = R^ + T 1 r 

s s 


and begin the manipulations with 


3 *s _ -*T 8T 

9^ k r 3 * k 


k = p, y, r . 


*See References 8 and 14 (among others) . 

**The well known special case of the generalized angular momentum equation would 
apply here, that being the case wherein the linear moment relative to the MRP 
vanishes identically. 
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Clearly, the manipulations can proceed no further until the matrix T is definitely 
specified by a prescribed Euler sequence. Fortunately, as extensive (yet simple) 
manipulations of the author have shown, the form (4-57) of the moment equation is the 
same no matter which of the 3! = 6 distinct Euler sequences below is selected to define 
the matrix T (see equations (3-0.1), (3-0.2), (3-0.3)). 

^y'* (3) [ V(2) ^(1) (a 1,2,3 sequence through «P p , and >P y ) 

[ *^P ] ( 2) tl V(3) ^r ] (1) (a 1,3,2 sec l uence through <P r , f y , and >P p ) 


T = 


[^pl( 2 ) (a 2,1,3 se q uence through * , «P p , and <P y ) 

[ ^ r ] ( i) tf ] (3) [^ 1 ( 2 ) (a 2,3,1 sequence through ^ , and <P p ) 

[f r ]( 1) t' / ’ p ]( 2 ) ^y ] (3) (a 3,2,1 sec l u ence through »p , >p p , and v? r ) 

[^p] ( 2) [ V(1) [ ^y ] (3) (a 3,1,2 sec l uence through «p , «P p , and >p p ) 


Selecting the 2, 3, 1 sequence (merely for the sake of being definite) so that 


r 


" r " 


I 



s 


s 

-> 




-> 

3 

= T 

J 

s 

[ V(1) [ V(3) [ V(2) 

J s 

_k_ 


-> 

_K _ 

s 


-K . 

s 


it can be shown with little effort that 


9 R 

£ 

3 


r T T 


0 0 1 . 

0 0 0 

- 1 . 0 0 


9 R 

< 

9 


?T T l V<» 


LO 


o -1. o' 
1. o o 
0 0 


[ V(2) 
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8 R s 
9 <P r 


' ^ 'Vd) 


0 0 0 

0 0 - 1 . 

0 1 . 0 


[ ^y ] (3) [ V(2) 


Further scratchwork (not repeated here) has shown that summing corresponding mem- 
bers of the three equations whose forms are those assumed by equation (4-1) when 


q = , v , and <p , followed by some rearrangement, gives 

3? y ^ 


1 + <V>B + < J s>B 


/ 


(m) 


r x T R s dm - f r X T (d ? GR % R T T Y H) S ) 

(m) 


J ' rx T Ic 

Cm) ' 


(SUN) 

GRAVITY 


,s) - / ?II ( dF 

(m) 


(MOON) 

GRAVITY 




J"r x T <d F AERO s ) - y? x T (d F SOLAR s ) |=0. 
(A) (A) 


(4- 54) 


In equation (4-54), the subscript B indicates the B -resolution, that is, 


(J ) n = B -RESOLUTION OF J = T J = T 
s B s s 


'O' 

1. 


(K s’>B 


[ ^(1) tV(3) K s’ 


[ ^r ] (l) [ ^y ] ( 3) 


= B -RESOLUTION OF K ' 

s 


where 


—y _ 

ri ' 


—y 

ri i 

s 


s 

—y 


-> 

J ' 
s 

li 

i—i 

to 

j 

s 



~y 

-K lJ 


Lk J 
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Recalling the definition of T = T (S+B), it should be obvious that 


J r x T d(§ 

(m) 


■ (EARTH) 
GRAVITY 


* 2 ^ = fvx d Fgg = Mg B [see (A- 28) and (A- 29)] 


(m) 


= TORQUE DUE TO EARTH'S 
GRAVITY FIELD 


/ 

(m) 


r x T d^F (SUN) ) 
r x i u \ r GRAVITY ,s/ 


M^ UN) = TORQUE DUE TO SUN'S GRAVITY FIELD 


/ 

(m) 


rxT d(l^ (MOON) ) 
r x 1 a \ r GRAVITY ,S/ 


= M^g 00N) = TORQUE DUE TO MOON'S GRAVITY FIELD 


Jv x T (d F aerq g ) = M aerq = AERODYNAMIC TORQUE (APPENDIX B) 


(A) 


J' v x T (d F s0LAR s > = ^SOLAR = SOLAR RADIATION TORQUE (APPENDIX C) 


(A) 


Treating both Sun and Moon as point masses results in the following approximations 
to Mgg^^ and these being found by developments similar to that of 

equation (A- 7). 


( SUN) 


3 y 


gB 


as 


R 


U SUN X ( 


Du, 


SUN- 


SUN 


M (MOON) 3 U M 

^ -,3 


gB 


R 


"MOON x ( D"i 


MOON 


) 


MOON 


(4-55) 


(4-56) 


In equation (4-55) 


R 


SUN 


= R 


SUN 


R 


SUN 


= T (R 


CM 


. £(SUN) } 


u 


SUN 


r sun /r sun 


while in equation (4-56) 
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r moon r moon 


p = T (R CM - R (MOON >t 

R MOON T (R s s ' 


u, 


MOON 


“ r moon /r moon 


The symbol LI denotes the inertia matrix of the entire system (in its instantaneous 
deformed configuration) referred to the B -frame which has origin at the instantaneous 
CM of the entire system. It follows from equations (4-54), (4-55), and (4-56) that 


/ r x TR dm = M „ 
s gB 

Cm) 


+ £ (SUN) 

gB 


+ M 


( MOON) 
gB 


+ M + M 

AERO m SOLAR 


(4-57) 


A strong advocate of the direct use of either the principle of angular momentum or 
equation (4-57) might insist that equation (4-57) be invoked at the outset in complete 
disregard of the manipulations leading to it. The author of this report would then 
hasten to point out that some of the readers of this report are likely to be the 
uninitiated who will find it instructive to base the development of the equations of 
motion upon the principle of virtual work. 

Alternatives to equation (4-1) in the approach to the moment equation are the 
relations (4-58) and (4-59), both of which are established in Reference 1 (see also 
Reference 15 on quasi-coordinates) . 


/ 


-"CM r 4 -» 

(KE) dm + (TRg ) x / r dm = M £XT , 

(m) (m) 


(4-58) 


f L^ (KE) dm + V x j V v (KE) dm = M EXT . 
(m) (m) 


(4-59) 


The operators and V v are defined by 


-( 8 ), 


L ( ) = 1-^1 V ( ) + w. „ x V ( ) = V ( ) + q T V ( ) 

w ' dt 'B w B to \dt^ B to v to 
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where 


V 

v 


9V, 


9V, 


-1 




-> 

+ k 


9 

9V 3 


9 

9V 3 


J (B) 


V 

0) 


and 


3 

3 co ^ 
3 

9u> 2 




9 

9 a) 3 


9 

9o) 3 


(B) 


V = 


V, 


V 


= V.i + V 9 j + V,k = TR 


CM 


(B) 


In both equation (4-58) and (4-59), the symbol KE denotes the kinetic energy per 
unit mass. However, in equation (4-58), it is necessary to write 


~ -+t> -» \ ->rM t -* 

KE = (1/2) R 1 R = 1/2 <R + T 1 (r + 


r) 


T f 4 


R 


CM 


+ T (r + fi r) 


while in equation (4-59), one must write 

KE = (1/2) jv + r + r 

— ^ 

The vector M is the B -resolution of the resultant external moment about the 
hx i 

system CM , and if comprised of the moments of only those external forces accounted 
for in this paper, it is identically the right member of equation (4-57). Complete 
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equivalence of the left members of equations (4-57), (4-58), and (4-59) is also estab- 
lished in Reference 1. 

-> 

Equation (4-57) with its right member replaced by is the same as equation 

(17) of Chapter 3 in Reference 8 wherein the symbol a denotes the acceleration dis- 
tribution instead of TR . 

s 

On writing 


TR 


s 


TR 


CM 


+ r + 2 W 


r + (ft‘ 


+ ' 


2 

) r 


and recognizing the relations 


/ 


r dm 


0 (by definition of the CM) 


i 


r x si' 1 ’ r dm = []u 


B 


J 


r x Si r dm = 




= Dw 


B 


J' { r* x Sp v + r x ?} dm = []t 


B 


(m) 


f 

(m) 


,T 4 


-□ 


2 /rxfi A rdm=]|aTg + (JgX / r x r dm = QuTg + ^ 


/ 


fi' 


(m) 


/ 


(r x r) dm , 


(m) 


use being made of the vector identities ax [b x (b x a)] = -b x [ax (a x b)] and 
a x (0 x c) = (a*c)b - (a*b)c in establishing the third and fifth, one can show that 
equation (4-57) may be put in the form 


D* B + tk B + ^ D“ B + f r x r dm + fi' 

(m) 


/ 

(m) 


r x r dm = Mg XT . (4-60) 
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Before dealing with the two integrals in equation (4-60) in detail, their integrands 
will be transformed via the substitution r = r - r^,^. Thus, the moment equation will 
be rewritten, temporarily, as 


O^B + Q^B + qT + f 

(m) 


(r - r CM ) x r dm + f(r - r CjVJ ) 


x r dm = M 


EXT 


(4-61) 


where, obviously, the condition jv dm = (f has been imposed. 

(m) 


To express the two integrals in equation (4-61) in terms of the system coordi- 
nates and known system parameters, the system will be subdivided in the manner 
indicated by the subscripts on the integral signs in equation (4-16) and substitutions 
made from equations (4-21) through (4-43). Manipulations leading to the results that 
follow will not be repeated here. 

It obviously follows from equation (4-21) that 


(m ) 
v o' 


fir - r CM ) x ? dm + Q 1 f (? - r CM ) x r dm = (f . 


(4-62) 


“V 


On substituting from equations (4-22) through (4-24) and again drawing from 
References 1 and 8, the author has arrived at the approximation 


/ 


r x r dm % - 


(m f ) 


/ 


r x ? (p r • n) dA % -M 


T 5 


(4-63) 


( a e 1 


~ ► 

the symbol M T denoting the moment of the thrust about the system CM and defined by 


/ 


-> r* 


= / rxr(pr-n)dA 


(A e ) 


•«/ 


r x (p 


P 0 ) n dA 


(4-64) 
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The approximations (4-63) are consequences of the complete neglect of the volume 

r x gj (p r) dV and the rightmost surface integral in equation (4-64), 

(V f ) 

the volume integral being encountered in the manipulations pertinent to J r x P dm. 
Also made in this report is the approximation (mp 


f r x r dm % 0 . (4-65) 

(mp 

In each of the vector products appearing in equations (4-63) through (4-65), the 

— y — >■ — ► 

factor on the left has been written r instead of its equivalent r - r^^ to better 
emphasize the fact that the moment is about the system CM. 

By referring to equations (4-25) through (4-27) and to the definitions and 
relations of sections 2 and 3 pertaining to swivel engine i, one can show (as in 
Reference 1) that 


/ 


(r - 


r CM ) x r dm + Q.' 


( m Ei> 


/ 


(r - r CM ) x r dm = - m E . q £i x Ag- 


( m Ei ) 


,T T Ei 


m Ei £ Ei W B X ^ q Ei X A Ei^ + T Ei 1 T Ei ^Ei + ^ W B + w eP 
X T Ei ^ T Ei "Ei > 1 = ••• NSE • 


(4-66) 


where 


q Ei r Ei r CM ^Ei A Ei * 

1L 

A relation similar to equation (4-66) holds for the i in rotor, as consideration 
of equation (4-35), equation (4-36) and the definitions and relations of sections 2 and 
3 concerning rotor i will show that (Reference 1 again) 


/ 


(r - 


r CM ) x r dm + $■' 


f x ? dm ~ T Rj 1 T 


( m Ri ) 


(m Ri ) 


+ <« B + “rP x tL I Ri T 


Ri w Ri 


Ri " Ri Ri 
i = 1 , ... NR 
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(4-67) 



From the content of the paragraph containing equations (4-28) through (4-30), 
there follows 


/ 


(r - r CM ) x r dm + fi' 


/ 


(r - r CM > x r dm = m pj «r pi - r CM ) 


<m pi ) 


(m p .) 


X (5 Pi A Pi } + U B x [(? Pi “ ? CM } x (5 Pi A Pi^ } , i = 1, ... NP . 

(4-68) 

Substituting from equations (4-31) through (4-34) and manipulating in accord- 
ance with the relevant definitions and relations of sections 2 and 3 leads to 


/(? 

(m.) 


r CM ) x r dm + ft 


J '^ r x r dm T i 1 T i °i' + + “i*^ 


(mp 


x T? I 1 T. up + m. (?. - r CM ) x T? [(fi? + ^ ) (I (0) + t (l) n (l) ) 
+ 2 ti? Y (l) n (l) + ^ (1) n (l) ] + w B x {m. (?. - r CM ) 


x ( r. (0) + ; (i >) + ;w ]} 


N. 


N. 


? i T tM / ? i x dm + n J XX / ** 


^ X dm 


j=l I (m.) 


K=1 (m.) 


N. 

l 


2 f * i X ( ^i X dm + 2 n jX^K f ^j (l) x ( “i x 


dm 


(m.) 


K=1 (m.) 


N. 

l 


+ w B x < T 


n: 


/ ? i x ? j 


(i) 


j=l L (m.) 


N. 

l 

dm + nj f ^ x dm 

K=1 (m.) 


i = 1, ... NA 


(4-69) 
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where 



and the symbol I 1 denotes the inertia matrix of flexible appendage i, in its instan- 
taneous deformed state, referred to the axes x.y.z. defined in section 2. 

ri 1 

For convenience in writing the moment equation below, the integrals in equation 

-*i ->2 i i 

(4-69) will be expressed more compactly by introducing the symbols C ., ., 

J J J 

and •$ ^ defined by 



/ 


r. 

l 


->(i) 

X 


dm 


(4-69.1) 


C V /V”* 5 K <1>dm 

(mp 

(i = 1, 

J { r] = f ^(r.)^(^ (i) ) dm 

(m.) 


(4-69.2) 

NA, j = 1, ... N., K = 1, ... N.) 

(4-69.3) 




(4-69.4) 


(m.) 


j -*• t 

where the operator P is such that when applied to the vector A = [A^, Ag> A^] , 

the result is the skew symmetric matrix 


Aa) = 


0 

A, 

-A„ 


-A, 

0 

A. 


(4-69.5) 
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The integrals in equations (4-69.1) through (4-69.4) also appear among the time 
independent integrals of Appendix F which arise in deriving the bending equations. 
Definitions (4-69.3) and (4-69.4) permit one to write 


f x (Sfj x £ (i) ) dm = -J r ] w. , f £ (i) x (1 x ^ k (i) ) dm = -J. ^ . 


(m.) 


(m.) 


The contribution of the i in SDOF CMG to the two integrals in equation (4-61) 
is, in view of equations (4-37) through (4-40) and (3-9) through (3-25), 


f (r “ ?cM ) x ? dm 


T / =i 

+ ^ J (r r CM ) 


x r dm 


<m Gi’ +m gi'> 


(m^.'+m .') 
Gi gi 


= (1 xx + >£> [ =Gi’ + “B x “Gi’> + 't t<“B + “Gi'> x “gi’ 1 ' 


(4-70) 


The right member of equation (4-70) is a simplification of the expression (4-71), 

< T Gi ‘ Gi T Gi + T gi lg ‘ T gi> “Gi + <“B + “Gi’> X (T Gi * Gi T Gi + T gi ^ T gi> “Gi' 

+ (!?„+« .') x T T . I gl T . u + 2 U .* x T T . I gl T . U ' 

B gi gx gi gi gi gi gi Gi 


+ [Tr (I gl )] (w Gi * x <? .') , 


(4-71) 


the first three terms of which one might have anticipated after inspection of equation 

(4-67) and certain terms of (4-66) and (4-69). By Tr (I gl ) is meant the trace of 

I gl , that is, Tr(I gl ) = I gl + I gl + I gl . 

’ v ' xx yy zz 

In this paragraph , as in that segment of section 3 pertaining to a 2 DOF CMG 
and in that paragraph of section 4 containing equations (4-41) through (4-43), an 
additional subscript i (or superscript i as the case may be) will not be attached to 
any symbol relating to a 2 DOF CMG. Manipulating in accordance with the definitions 
and relations of the aforementioned segments of this paper leads to the contribution 
of a 2 DOF CMG to the two integrals in equation (4-61), that being 
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f 


(r - r CM ) x r dm + 


^ nl OG +m IG +m g ) 


/ 


(r - r CM ) x r dm 


<m OG +l,, IG +m K ) 


T OG+IG+g -5- , , ^ T OG+IG+g -*• , 

" *B W OG + + W 0G ^ X l B w OG 


+ lo g ^ w Ig' + “iG* X ^OG^ + ^“B + W IG ? ^ 


B 


x I 


I G+g 


(I)./,' + 2 a> T ’ x I 


I G+g 


B IG IG A B OG 


- [Tr (I IG+g )] (u IG ' x u 0G ') + (u fi + «;) x I B g u • 


+ 2 “V x ‘I <“og' + “ig'> 


- [Tr (I s )] [JT ’ X (JT og ’ + u IG ')] 


(4-72) 


The factors Tr(I^), E, = IG + g, g, in equation (4-72) could be replaced by Tr(Ig), 
5 = IG + g, g, since the trace is invariant under a similarity transformation. See 
equations (3-28.1). 


On assembling the several contributions to the integrals in equation (4-61), as 
given by equations (4-62), (4-63), (4-66), (4-67), (4-68), (4-69), (4-70), and (4-72), 
one can recast the moment equation in the somewhat more familiar form 


O^R 


•B ' LJ"B + 3 B x 0"B = “l + “gB + MgB° 0N> + * 5 gI UN) + “AERO + “SOLAR 


NSE 


+ m Ei £ Ei * q Ei X A Ei + W B X ^ q Ei X A Ei^ 


i=l 


NSE 


X] ^ T Ei 1 T Ei w Ei + ^B + w Ei ^ X T Ei 1 T Ei w Ei * 
i=l 

(4-73) 

(continued on next page) 
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NR 


Z < T Ri I Ri T Bi » Ri ' + (« B + « R1 ') X T l t I Ri T Ri J Ri ’} 


i=l 


NP 

m Pi {(r Pi " r CIVl ) x ( ^Pi A Pp + w b x [(r Pi 
i=l 

x (5 W A pi )]} 

NA 

L {? i T ^ *i “i’ + + “i'> X ^ ? i “i’> 

i=l 


r CM ) 


NA 

Z) < m i (r. - r CM ) X T? [(fi£ + flj ) (r. (0) + Y (i) ? (i) ) 

i=l 


+ 2 v ^ n ( - 1 ' ) + n^^]} 


NA 


B 


x £ {mj (r. - ? CM ) x (r. (0) + f (l) n (l) ) 


i=l 


+ * (i) ^ (i) ]} 


NA N i 


N. 

1 


i=l j=l 


K=1 


2 ru J } u. 

3 rj i 


N. 

1 


- 2 nj £ n^i* J) 


“b x 




ii 


NA 

N. 

( i 

N. 

r i 


E 

tez 

pi c oj + "i E n R 

c 1k] 

i=l 

' i=i 

K=1 



(4-73) 

(continued on next page) 
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NSDOF 


E 

i=l 


« + ‘S > «Gi + "B X ^gP + l t <“B + ”gP 


X (I)' .} 

gl 


N2DOF 

I 

i=l 


E a° G+IG+ ® i bc + ( » B ♦ s bG ) x i GG+IG+ e » bo 


. T I G . — > . ~> # ~> 4 x — ^ . 

+ (“In + x + (“n + “in) 


IG 


IG ~ OG" 


B IG ; 


x I^S ? w ♦ 2 3J 0 x ;< 0G 

- [7r (I IG+g )] (^ G X ^ G ) + (<^ B + 3* ) 

x l| ^ + 2 3* ff x I§ (JSJjq + 5P iq ) 

- [7r (I g )] [w* x (uq C + w' )]} . (4-73) 

g «> (Cone.) 

T Ei -> T Pi 

One expecting to see terms such as -T* . I T„. <*>„ , -T„ . I T 0 . u- etc., and 
T Ei _> t Ri 1 111 ^ Ki B 

-ojg x Tgj I Tg. JTg , -ojg x Tgj I T R j oTg , etc., in the right member of equation 

(4-73) should not be alarmed since their negatives are embedded in the left member 

of equation (4-73). Recall that [] is the inertia matrix of the entire system referred 

to the B-frame defined in Section 2. An expression forQis given in Appendix E. 

The primed vectors associated with the CMG’s in equation (4-73) have the 
B-resolution as do the vectors ojg. , Wgj, and to! (see Section 3). Their introduction 
made it possible to write the moment equation more compactly, not to mention impart- 
ing a certain "consistency of appearance" in the presence of the other primed vectors. 
The explicit dependence of the CMG terms upon the gimbal angles and their time 
derivatives and the gyro element spin rates is shown by the following: 


(I S + ■£> <“Gi + “B x “gP = (I S + '£> 


*01 T L 


+ 6 Gi W B X l T Gio 
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l fz ( “B + “bi } X “ gi ' 


-IS 1 w 
zz 




gi l* 5 Gi ^Gio 


cos 6 g . 


sin 6 


Gi 


" “ B X l T Gio 


OG+IG+g 
a B OG 





L 

J 



r (I g - 

ZZ 

I g ) 

yy 

sin 

6 og sin 

6 ig cos 6 ig 

^ OG T BGB 

(I g - 
yy 

i g ) 

Z Z 

cos 

6 oa si" 

6j G cos 6 iq 


lOG + 
ZZ 

i IG + 

ZZ 

i« 

yy 

■ 2 . 
sm S IG 

+ 1^ cos^ 
z z 


0 

-sin 6 
cos 6 


6 ig j 


Gi 
GiJ 
(4-73.2) 


(4-73.3) 


, \ T OG+IG+g ->t _ -► v 

(u fi + “0G ) x X B w OG W B 


| ^ OG T 


T 

BGB 


(1* - 
zz 

I g ) 
yy 

sin 

6 og 

sin 6 j G cos 6 

(l g - 

yy 

i g ) 

Z Z 

cos 

6 og 

sin 6 j G cos 6 

lOG + 
ZZ 

I IG + 

ZZ 

i* 

yy 

sin^ 

6 IG + l lz C0S 


IG 

IG 

2 


^IG"^ 



ra g - 

zz 

I g ) 
yy 

cos <s 0G 

sin 

6 IG cos 

6 IG 

. r 2 T T 

0 OG BGB 

(I* _ 
zz 

i g ) 
yy 

sin 6 oq 

sin 

6 IG cos 

6 ig 




0 



. 


l B g (w IG + W IG X W OG' ) 


(I IG + I g ) 6 

v XX XX 7 


ig t bgb 


cos 6 


OG 


sin 6 


OG 


(4-73.4) 


(4-73.5) 


(w B + W IG ) X ! B S “IG 


(I IG + l g ) <S T „ 
v xx xx ; IG 


T 

,“B BGB 


cos <5 


OG 


sin 6 


OG 


(4-73.6) 


53 


sto S 0G I«i2 * ‘yy < 2i IG> 


+ 2 I s sin 6 , n ] 
z z IG J 


{2wJ g x I B g w oG [ 7r(I g )J (wj G x ^q Q ) } 5 qG <5 ig T bgb cos 6 qg [I xx + I yy cos ( 2<5 IG ) 

* 2 lf s sin 2 t IG ] 


“fy ' 'lz> si " < 25 IG> 


(4-73.7) 


-sin 6q G cos 6j G 

( W B + x *1 w g = r yy w g | W B X T BGB COS 6 OG COS 6 IG 

sin 6 T „ 


(4-73.8) 


(2^ x l| (^ n + u ’ r _) - [Tr(I® ) ] [< x ( W ^ G + o)» G )]} 


g B OG IG 


cos 6, 


“ w g *yy t bgb roc cos 6 ig sin 6 og 


-sin 6 j G sin 6 GG 

sin <5 j G cos 6 qg 

-cos 6 T/ _, 
lu 

(4-73.9) 


The relations (4-73.1) and (4-73.2) hold under the assumption that I Gl and I^ 1 are 

diagonal and 1^ = I^ 1 , while equations (4-73.3) through (4-73.9) are valid if I GG , 

I IG , and I s are diagonal with I 0G = I OG , I IG = I IG , and I g = l g 
’ & xx yy ’ yy zz’ z z xx 


In view of equation (4-53), it seems that the selection of the system CM as the 
MRP, rather than the origin of the xyz axes (the B-frame), in deriving the moment 
equation, (4-73), is an inconsistency. A more practical and wiser choice of MRP could 
be the origin of the B-frame, a point to be discussed subsequent to writing equation 
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(4-73.14) below. Equation (4-73.14) is the form of the moment equation when moments 
are referenced to the origin of the B -frame. The author’s derivation of (4-73.14) 
begins by invoking (4-57) which is here rewritten more concisely as 


y? x TR s dm =^eXT (4-73.10) 

J^CM 

EXT ^ en °tes the resultant moment of the external forces about the 
system CM (and has the B resolution). Via the relations 


r = r - r 


CM 


r CM = f r dm (see Appendix D) 

(m) 


-*■ T~ . T T ~ 

TR = TR + r + 2£2 X r + (fi + fT ) r 

s s 


one may pass from equation (4-73.10) to (4-73.11) 


f r x (F + 2£1 T r + (si T + q T ) r} dm = it + * 

(m) 


CM X J' ir + 2Q r 
(m) 


+ (^ + ) r} dm 


(4-73.11) 


With denoting the B resolution of the sum of the external forces, it follows 

from equation (4-73.11) and either equation (4-44) or (4-45) that 


/ zz ^ t ^3 . t T ^ u CM ~ 

r x {r + 2fi r + (ft + ft ) r} dm = JL EX -p + r CM x ^ £XT - tn r CM x T R g 


(m) 


(4-73.12) 


But 


£CM ~ 2 Jjf 

‘M* EXT r CM X -^EXT *=%XT 


= sum of external moments about origin of xyz 
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r x si' 1 ' r dm = [k B > (0 = inertia matrix of entire system referred to 

xyz) 


B x Qw B = Q( 


r x Q r dm = x = Q 


2 f r x r dm = QuL + Q' 


x r dm 


so that equation (4-73.12) may be written 


D^r + Q 


B u B 


Wn + n T |_Ja) R + m T R = J/L 


s ^EXT " f r 

(m) 


x r dm - Q ' 


Dealing with the integrals in the right member of equation (4-73.13) in the manner of 
arriving at equations (4-62) through (4-69) and (4-70) through (4-72), there follows 


□4 +0 


UJn + Q (J0 T 


+ m T R g 

XT ^ e/ ^ > T 


+ {other terms attributed to the integrals in (4-73.13)} 


The vector is the resultant engine thrust moment about the origin of the B-frame, 
while the terms alluded to within braces include all the sums appearing in the right 

member of equation (4-73) modified, where necessary, by suppressing the symbol 

— > 

?civj. The vector q^. in the sum pertinent to the swivel engines should be replaced 
by *^Ei “ *Ei ’ £ Ei A Ei ( reca11 that ^gi = ^gi “ *CM " A Ei A eP" 

If [] is known, then one can compute G in accordance with the relation 

G = D+m^(r~ CM )£(-? CM ) . (4-73.15) 

However, the author does not recommend the use of equation (4-73.15) even if equa- 
tion (4-73.14) is chosen to be the moment equation. With little effort, equation (E-4) 
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of Appendix E can be transformed into one for determining G directly. Such a 
transformation* is effected by merely replacing the symbol I f with T f , the symbol T f 
denoting the inertia matrix of the fluid referred to the B-frame, and replacing each 
symbol by the symbol where q^ is defined by 


q 5 q S + r CM 


K = 0, i, Ei, Ri, Pi, Gi, gi, OG , IG , g 


(4-73.16) 


the qj. (5 = 0, l, Ei, Ri, Pi, Gi, gi, OG, IG, g) being exactly as defined in appendix E. 
* _► % 

From the definitions of the q and the q , it should be obvious that the expressions 
for the q^ are devoid of dependence upon r^, and hence, that the resulting expres- 
sion for Q will be neither explicitly nor implicitly dependent upon . A differentia- 


I I py 

tion with respect to time would then yield an expression for M which is free of both 
rcjyj and r^^. Computation wise, the use of such expressions for [] and Q would be 
far more expedient than the direct application of equations (E-4) and (E-5) followed 
by an application of equation (4-73.15) to find Q and [] In fact, computation of [] 
and 0 would proceed more efficiently by first computing Q and Q and then [] and \\ 
from 

Q = D" m ^(r c]V ]) 


D=U- m £(? cm ) £(-? cm )- 


m ^” r CM^ " m 


From a programming standpoint, it should be evident that the choice of the 
origin of the B-frame as MRP offers certain computational advantages. In that regard, 
still another point favoring that choice should be made. Consider recasting the 
system equations of motion in the form 


AX = F (t, X, X) 

as a first step in rendering them amenable to some numerical integration scheme. (As 
one submatrix of the partitioned column matrix X, the author has in mind the 3x1 
column tjTg = [i^, ijj 2 , ijjg] whose elements are the second time derivatives of the 

*The reader should realize that here the author is not changing the appearance of 
(E-4) via the substitution of q - ? CM for q ? , $ = 0, i, Ei, ..., g. If such a 

substitution were made, after replacing I by I , the resulting right member of (E-4) 
would define neither [] nor Q ; however, if I f is left intact, the substitution would 
result in an expression which still defines Q, not []. It should be clear that the result 
of making the "replacements” alluded to above is the same as that obtainable by invoking 
the "generalized transfer theorem" expressed by equation (E-3). 
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quasi-coordinates tp. = J id. dt, i = 1, 2, 3, which is to say, obviously, that = Wg . 

Also, the author has here supposed that equation (4-74), not (4-76), governs the 
motion of the point mass m^,..) With the combination of equations (4-53) and (4-73) 
designated as system translational and rotational equations, it will be found that the 
mass matrix A is non- symmetric! But, that the matrix A corresponding to the com- 
bination of equations (4-53) and (4-73.14) is a symmetric matrix. The time consumed 
and storage space required for constructing and inverting the symmetric matrix should 
be less than that for the non- symmetric matrix. The author would definitely prefer 
equation (4-73.14) to equation (4-73). 


Clearly, the moment equation must be accompanied by the equations giving the 
Euler angle rates, , K = P, Y, r, in terms of the components of uTg and the ^ , 

K = P, Y, r. The structure of these equations will, of course, depend upon the 
"Euler sequence" of rotations selected to pass from the S -orientation to the B- 
orientation. For a 2, 3, 1 sequence through the angles ^p, Py , respectively, 

the ^ K , K = P, Y, r, are determined by equation (3-0.4) except in those cases where 
v? y is an odd integral multiple of tt/ 2, thereby making the matrix in equation (3-0.4) 
singular. Attention will not now be diverted to methods for circumventing the 
numerical difficulty encountered at such singularities. 


With the coordinate 5 . assuming the role of q in equation (4-1), the result of 

some manipulation is the equation governing the motion of the point mass m . , 

Pi 

m Pi ^Pi + C Pi ? Pi + K Pi ? Pi " m Pi l w B x A Pi 1 ? Pi + m Pi ^Pi {T (R s “ ^Pi ) 


+ (si T + 0, T ) r p .} = 0 


i = 1, ... NP . 


(4-74) 


The only symbol in equation (4-74) not heretofore defined is A^. which is given by 


A . 
Pi 


3) C ' 6 pi ] ( 2) (R pi’ X 


l . , 6 .) + u 

pi pi ! 


- (SUN) _ £ pi 
(SUN) 

I AV c 


- R P ‘| 3 

s 


R (SUN) 


^ (MOON) 

- R P 1 

\ v 

£ (MOON) 

n s 

-U i | 

s 

s 

s I 

IR < SUN) I 3 

y M j 

(MOON) 

' s 

- R pi l 3 

(MOON) ,3 1 

1 s 1 1 

i 


i 

s 1 

' s ' J 


(4-75) 
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The definition (4-75) must be complemented by those following: 


R P1 = [X p \ Y*\ zn* = R + T A (r . + 5 . A .) 
s L s ’ s s s v pi ^pi pr 


,pi vP 1 7.P*i T - ^ . mT 


= position vector, referred to the S- frame, of m . 



“pi = tan_1 < Y f /X s i) 


right ascension of m . (0 < a . < 2ir) 
& pi v pi 


6 . = sin 1 (Z pi /R .) 
pi s pi 7 


declination of m^ (-ir/2 ^ <5^. < ir/2) 


A . = 

Pi 


°pi “ “p " w^(t~t o ) = east longitude ofmj(0<A.< 2ir) 


(If the expression for Ap. is negative then A . is to be replaced by its positive 
equivalent modulo 2 tt . ) 


A (R A ., 8 .) = A evaluated at point with spherical coordinates 

o Xr* 1, ir r ' 1 o 

(R ., A ., 8 .) referred to the E -frame (Section 2) 

r* Jr 

(See equations (A- 2) through (A- 5) of Appendix A and 

that part of Section 3 on the unit vectors u D , u. and u_ . 

K A o ) 

It will be observed that in writing equation (4-74), the vector r . was regarded as a 
constant vector as implied by equation (4-29). 

To allow for the possibility that m . is a "trim" mass installed for attitude con- 
trol, one could add to the right member of equation (4-74) the force F^. /T., the 
scalar function of time F . being determined by some control law. 

2 

Introducing the symbols to^ and , satisfying m^. and 2 m^ r . co . 

= C and denoting, respectively, the undamped natural frequency and critical 
damping ratio associated with m^., permits one to rewrite equation (4-74) as 
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(4-76) 


r . + 2 r . oj . £ . + (oj . 
^pi h pi pi ^pi v pi 


l“B X ? pi! 2> «pi + 'pi {T (S s - V 


rp rp " 

+ (n + n ) r pi > = o 


i = 1, ... NP 


Despite the title of this report, the influence of gravity on the motion of swivel 
engine i will be accounted for only in part , while that of the aerodynamic force and 
solar radiation pressure will be neglected entirely. In approximating the effect of the 
gravity field, the strength of the field is assumed constant throughout the region 
occupied by m E - and equal to that at its CM. The resulting equation corresponding 
to the engine deflection 3 . is, before further simplification 




9 3 


Pi 


^ T Ei 1 T Ei + w Ei + W B x w eP + ^ W B + w eP X T Ei 1 T Ei ^B + w eP 


3 A 


- m„. 


Ei 


Ei Ei 3 3 . 

Pi 


T < R s - A Ei> + 


(SI, 


<r Ei ' I Ei A Ei ) 


+ C pEi »pi + K pEi <B pi ' 6 PCi> 0 ■ i = 1. ••• NSE . (4-77) 


The companion equation, corresponding to 3yj> can be written immediately by merely 

replacing the subscript P in equation (4-77) with the subscript Y. (This is not to 

say that the final form of the equation satisfied by 3y. can be obtained from that for 

3pi by replacing the subscript P with Y. It is only at the stage of development 

indicated by equation (4-77) that such a change in subscript may be made.) In 

writing equation (4-77), the author has made use of the fact that 3R /3q = 3R /3q 
• s s 

and 8A Ei /83 pi 9A Ei /93 pi- The equation defining A £ . is quite similar to equation 
(4-75), that being 
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where 


o Ei _ r „Ei Y Ei Ei T _ 5 , t T ( ~t , 

R s = [X s ’ Y S ’ Z s ] " R s + T (r Ei " £ Ei A Ei } 


= position vector, referred to the S-frame, of the CM of m 


Ei 


R Ei= 1 * 5*1 


ctgj = tan 1 (Y R */X R *) = right ascension of CM of mg. (0 ^ < 2 tt) 

<$Ei - sin * (Z Rl /Rgj) = declination of CM of (- ir/2 6g^ ir/2) 


Afii = Og. - Op - « (t-t Q ) = east longitude* of CM of nig^ (0 ^ Ag. < 2-rr) 


A g (R Ei , 


Afii > 5gj) = Ag evaluated at the point with spherical coordinates 
(Rg., Ag., Sg.) referred to the E-frame. 


Obviously, one can proceed no further until the components of K„. = [A„., y_,., 

rp XjI fcl lil 

Vg.] are given explcitly in terms of 0 . and 0 Y j, and such expressions cannot be 
determined until the engine deflection sign convention, actuator arrangement, etc., 
have been decided. However, it is not difficult to deduce that 









Ei 




(4-79) 


whatever the structure of the scalar functions defining the direction cosines Ag., Vgj, 
and Vg.. The symbol Wgj, denoting the angular velocity of swivel engine i relative 
to the B -frame, was first introduced in Section 3 where the author failed to point out 
that the Xg. axis of the Ei-frame (which is fixed relative to swivel engine i) is 
directed as the thrust vector and hence as A_.. 

iil 

*It should be evident that if the expression for Ag^ is negative, it is to be replaced 
by its positive equivalent (modulo 2ir). A similar remark holds for any symbol denot- 
ing the longitude of a point. 
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The only special case of equation (4-77) to be considered in this paper is that 
which applies to an engine whose thrust vector, Fg., makes an acute angle with the 
x-axis, whose pitch actuator causes the xz-projection of the thrust vector to rotate 
about an axis parallel to the y-axis, and whose yaw actuator causes the xy-projection 
of the thrust vector to rotate about an axis parallel to the z-axis. For such an 
arrangement, it is easily shown (as in Reference 1) that Ag^ is given by 


L Ei 


Ei ‘ 


' 1 . 
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Ei 

K Ei 

tan s yi 

Ei _ 


-tan 3 . 


L P pi j 


where 


Kgi = (1 + tan 


yi 


+ tan 


Pi' 


^yi ^yi + Y yEi + e y E i 
^pi ~ ^pi + Y pEi + e pEi 


Y Ei = en &i ne can t an gle (0 < Ygj < W2) 


Y yEi 


engine cant in yaw (a known function of Yg. and the polar coordinates 

of the point of application of F„., the pole to which the polar coordinates 

h l 

are referred being the point defined below), - tt/ 2 < Yy E j < tt/ 2 


A Eio " the line of action of F Ei When e yi = e pi = e yEi = e pEi = 0 


Pg. - the point at which line Ag. Q intersects the plane y = 0 

(This point need not be a point of the x-axis though such was the 
case for the early Saturn vehicles.) 

P Ei " the Projection of Pg.. upon the plane containing the point of application 
of Fgj and parallel to the yz-plane. 
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£yg. - engine misalignment in yaw (- ir/2 << << ir 


/ 2 ) 


Ypg- - engine cant in pitch (also a known function of Yj^ and the polar 
coordinates of the point of application of F^.) , - tt/ 2 < Ypgj < tt/2 


£p E - - engine misalignment in pitch (- tt/ 2 << £p E - << tt/2) 


The expressions for an< t 3pj in Reference 1 included terms representing the con- 
tribution of bending, it being assumed there that the carrier vehicle experienced 3D 
bending at the swivel point. In this paper, however, the swivel point, supposed a 
point of the central carrier, suffers no displacement due to bending by virtue of the 
assumed rigidity of the carrier. 

Again, drawing from the author's previous work (Reference 1), one has (bend- 
ing or no bending) for this special case 


3 ^Ei 

8 3 . 

Pi 


- a E . (1 u Ei )/A £i 
U - v Ei 2 ) (1 - 7 m 2 ) 


3A Ei 

3B pi 


~ t - 2 - 
v Ei (i y Ei ; 


y Ei v Ei (1 ' y Ei )/X Ei 
^ y Ei ^ ^ " v Ei ^ X Ei 


3 ^Ei 




yi 


y Ei v Ei /X Ei 


y Ei v Ei 


1. 


3A Ei 

3 3 • 

yi 


- y Ei (1 - v Ei ) 

(1 - P Ei 2 ) (1 - v E1 2 >A Ei 


“ y Ei v Ei (1 ' V Ei )/X Ei 


63 


Ei , ~ 2 

1 ' u Ei 


y Ei v Ei + y Ei v Ei 


y Ei y Ei 

• ft 

( 2 "Ei ^Ei 


A Ei v Ei " A Ei v Ei 

^ , 1 ~ 2s 2 

a - v Ei > 

■ A Ei "Ei 

A Ei y Ei " y Ei A Ei 


*Ei “Ei ' u Ei "*Ei 


and, with no bending at the swivel point. 
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Now , "if" the angles (3 . and gyj are sufficiently small as to permit the approximations 
tan % . % (? ., tan 3 ! . % % ’ and further, if such products as ^In 5,2 are negligible 

]j± jjj. y x y 1 

when £ 1 + £ 2 - 2 ’ ? and h being any members of the following list, 
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> i 


y Ei’ y Ei’ y Ei’ V Ei’ v Ei’ v Ei’ 3 pi’ 3 pi’ 3 pi’ 3 yi’ 3 yi’ 3 yi’ “l’ w 2’ “3* 


Up °°3 9 


then considerable manipulation will show that equation (4-77) becomes 
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+ m Ei £ Ei [ *B (3 pi + Y pEi + £ pEi ) + + ^Ei “l x Ei w 2 ] 


+ C PEi B pi + K pEi (3 pi " B pCi ) ^ 0 


i = 1, ... NSE 


(4-80) 


where 


[X|\ Y* 1 , Z= i f= T (R s - A e .) 


Similar approximations made in the equation corresponding to 3 y . lead to 

(I zz + m Ei £ Ei ) ( 3 yi + U) 3 ) + *xz “l + X yz ( 3 pi + a, 2 ) 

+ m Ei ^ Ei [Xg (B yi + Y yEi + e yEi ) - Yg - x Ei w 3 + z £ . Uj] 

+ C yEi 3 yi + K yEi ^ 3 yi ~ 3 yCp ^ ® ’ i = 1» ... NSE . (4-81) 

Since a CMG is a "small" component, it is here supposed that the effect of 
gravity on the motion of a CMG gimbal is altogether insignificant. With that in mind 
and with q = 6 ^ . in equation (4-1), there follows 
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th 

as the equation governing the motion of the gimbal of the i SDOF CMG. 

Writing the equations satisfied by the outer gimbal angle 6 q G and inner gimbal 
angle Sj G of a 2DOF CMG requires considerably more manipulation than that required 
to arrive at equation (4-82). The equation of motion of the outer gimbal is 
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(4-83) 
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while that of the inner gimbal is 


(I xL G + I xx ) ^IG ' *yy w g ^OG cos 6 IG + (I fz " I fy ) ^OG sin 6 IG cos 6 IG 
+ [cOS 6 OG> Sin 6 OG’ 0] T BGB {(I xx + I L> “B 

+ “b x [t big a IG + t big ^ b ]> 


-tS 


r IG - T g 


x yy sin 6 IG sin 6 og ' 6 og Sin 6 OG [I xx ' *yy cos (26 ig ) 
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+ 2 1^ cos" 
z z 


6 ig ] 


”‘yy “g 008 4 IG • (I ?y " 4 OG sin (24 IG ) 


t bgb “b 


+ C IG 4 IG - •*«,„ = 0 • <4 - 84) 

In equations (4-83) and (4-84), as in certain other paragraphs, neither subscript i 
nor superscript i has been attached to any symbol relating to a 2DOF CMG. 

Appealing again to equation (4-1) with q = 0 . and manipulating in the light of 
the definitions and relations pertaining to flexible appendage i, it is found that the 
coordinate 6. must satisfy 
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(4-85) 


In equation (4-85), the symbols F . and M . denote, respectively, the force and 
moment exerted by the Earth's gravity field upon flexible appendage i, the reference 
point for M^ being the instantaneous CM of the appendage (in its deformed state) . 
The vectors F . and M . are given by equations (A- 30) and (A- 31) in Appendix A. 
Expressions for the differentials dF AERO and dF SOLAR are ava ^ able fr* om Appendices 
B and C . The subscript A. on the integral signs indicates (as one should expect) 
that the integration extends over the surface of appendage i. The surface integrals 
are, in general, amenable only to numerical methods. On observing the absence of 
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time derivatives of r . , it should be recalled that r. is a constant vector on the B - 


frame to imply that !*. = (F = ?. and, hence, that (d z /dt z ) r. = (ft 1 + ft" 1 " 2 ) 

IX SI 1 . 

Equation (E-l) of Appendix E provides an expression for the inertia matrix I 1 in 
terms of that of the undeformed appendage and the mode shape functions and general- 
ized bending displacement coordinates pertinent to the appendage. By equation (E-2), 
one can compute the matrix []* necessary to the computation of both F^. and M^.. 

(The terms reflecting the dependence of F^. on [J 1 can probably be safely neglected, 
however, not all terms showing the dependence of M . on [J 1 can be discarded.) 

til 

To find the y bending equation associated with flexible appendage i, one 
invokes equation (4-1) with q equal to the modal coordinate ry The result of much 
manipulation is 
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(4-86) 


Definitions of the symbols C^. , and^j^ will be found in Appendix F. 

Of the several terms in the right member of equation (4-86), the generalized 

/ n* \ 

force Q > } , attributed to the Earth's gravity field, is of particular interest. In fact. 
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the primary objective of the assignment out of which this paper grew was to deter- 
mine the effect of gravity on bending. The author has made the following approxima 


tion to Q , 
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(Continued) 
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(Concluded) 


By definitions given elsewhere (Sections 2 and 3 and Appendix A) , the product 
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which appears as a factor of the first term in the right member of equation (4-87) 

should be recognized as the i-resolution of the gravitational acceleration due to the 

Earth at the point with spherical coordinates (R A 6 .) referred to the E-frame, 

oi oi oi 

that point being the CM of flexible appendage i. The reader's attention is here 
called to that paragraph of appendix A which contains equation (A-30). The dominant 
term of the several terms comprising the expression 


T 

-► (\ \ ~~ T T - * 

m ,i. V- 1 / r r f T , r i x rr a t r \ 

i v i T i T LU oi J ( 3) L ~°oi J ( 2) A g (R oi’ oi’ °oi ; 


IS 


-M E m, T , 

R 2 *j u Roi 
oi 


(4-88) 


_ 5 

as equations (A-3), (A-4), and (A-5) clearly show. All terms through those in R . 
were not retained as the author first intended. Many were dropped because they 
could not be expressed as the sum of products of functions of time and integrals 
whose integrands are independent of time. Still others were discarded because of 
being either too unwieldly or of questionable importance. It is the author's contention 
that the approximation (4-87) is more than adequate for any practical application, and 
in fact, that most people would neglect all terms with the exception of expression 
(4-88) or its dominant part. Retention of only the terms belonging to (4-88) is 
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equivalent to ignoring the variation of the gravity field over the appendage and 
assuming its magnitude and direction at all points of the appendage to be the same as 
that at the appendage CM. Definitions of the D's with subscripts and superscripts 
in equation (4-87) are among those of the time independent integrals in Appendix F. 


The gravitational fields of the moon and Sun, ignoring their variation over 


th 


appendage i, make the following contributions to the j bending equation 
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In writing the mathematical definitions of the generalized forces q (SOLAR) and 
q(AERO), Qne can only ag ^ ar ag the i n ^ e g ra i expressions. ^ 
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the differentials dF^ ERO and dFgQ LA £ (here supposed having the B -resolution) being 
those defined in Appendices B and C. Analytical evaluation of the integrals in (4-91) 
and (4-92), which are both geometry dependent and time dependent, is, in general, 
impossible. 

With slight modifications, the equations of motion in the preceding paragraphs 
can be made applicable to the configuration of Figure 2 provided the boom is regarded 
as rigid. The necessity of the addition of a few terms to certain of those equations 
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is a consequence of allowing each solar panel an additional rotational degree of freedom, 
that being the freedom to rotate about an axis parallel to the axis of the boom to 
which it is attached. The symbol y?. will here denote the angle through which the 
solar panel designated "flexible appendage i" in Figure 2 rotates about an axis parallel 
to the boom axis, which axis is parallel to the y.-axis when 0. = y. - 0. Introduction 

* y .+ ^ ~ t 

of the coordinate jp. naturally subjects the expressions for w. and ft. of Section 
3 to change, the revised expressions being* 

^ i ~ ^ V(l) ^ ^ (2) 
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The rotation matrix T^ is still defined by Tft = Cf ^ T^, it being understood now that 

Tj = [T(B->i)] 0 _ ^_q» and the relations (3-6) and (3-8) continue to hold. Equation 
l i u 

(4-73), the moment equation, remains valid as it stands. As for the other equations 

of motion, the appearance of only (4-53), (4-85), and (4-86) need to be altered, 

while (4-80) and (4-81) should be deleted since there are no swiveled engines (the 

RCS thrusters are fixed relative to the boom). To account for the elastic restoring 

force -K <p. and viscous damping force -C both of which resist a change in 

p i 1 ^ i 

¥>-, one should add to equation (4-14) the expression (1/2) and to (4-15) 

i v i 1 

the expression (1/2) C the index i having the appropriate range. 

. IP . 1 

1 1 


*As evident in the expression for -A, it is here supposed that passage from the "null" 
orientation of the i-frame to its instantaneous orientation is effected by a 2, 1 
sequence through <p ^ and 0.. 
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After modification, equations (4-53), (4-85), and (4-86) read as (4-53), (4-85)', 
and (4-86)', respectively. 
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At least one more remark regarding equation (4-86)' is in order, that being 
that it can be made to apply to a flexible appendage which is denied rotational degrees 
of freedom (such as 0^ and y?p by merely setting u>. and to 0 and replacing T. by 

T.. When 6^ and are identically zero, the rotation matrix becomes the 3x3 
identity matrix and T. is the same as T\ . The allied equations (4-87) through (4-92) 
apply without change provided that is replaced by when 6. = 0 and <£. = 0. 

The system coordinate must satisfy 
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where the symbol denotes the control torque applied to induce the appropriate 
change in #>.. 

It is important to note that equations (4-85)' and (4-93) were written for an 
appendage similar to that in Figure 1, but having two rotational degrees of freedom, 
and on the assumption that the solar panels designated "flexible appendage i" and 
"flexible appendage i + 1" in Figure 2 may move independently of each other, and 
further, that these equations hold whether both the 9. motion and the p. motion are 
executed simultaneously or 6. is held constant while <*>. changes and is held 
constant while 0^ changes. Obvious simplifications are possible if either 6. or is 
constant while the other changes (in such a case, equation (4-85)' becomes (4-85)). 

If flexible appendages i and i + 1 of Figure 2 constitute a single unit and always move 
in unison (so that 9. + ^ = -0., = - ^) , then equations (4-85)' and (4-93) apply 

to the combination (which should be relabeled "flexible appendage i") provided the 
symbols with subscript i or superscript i in those equations are pertinent to the 
combination . 

Terms of (4-53)' and (4-73) which are impertinent to the configuration of Figure 
2, such as those attributed to swiveled engines (not present on this station), must be 
deleted as should an entire equation corresponding to a non-existent system coordinate. 

If a module, which may be regarded rigid but is joined to the central carrier 
by means of a flexible attachment, is not adequately modeled as a point mass with 
only one displacement d.o.f. relative to the carrier then, obviously, an equation such 
as (4-74) does not provide an adequate description of its motion. Permitting such a 
module to have three translational d.o.f's and three rotational d.o.f’s relative to the 
carrier would require that six equations be appended to the system equations of 
motion, and further, that still more terms be added to equation (4-53), the system 
translational equation, and to (4-73), the system rotational equation. Nomenclature 
pertinent to this type of component follows. 


Nrb = number of rigid components with 6 d.o.f. 's relative to carrier 

m RB = mass component with 6 d.o.f. 's relative to carrier 

i 


r RB = P osi tion, referred to B-frame, of the origin of x RR y RR z pR (RBi frame) 
i 


RB. y RB { RBj 
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x RB y RB z RB - a right-handed rectangular frame, arbitrarily but sensibly 
1 1 1 oriented, at rest relative to the B -frame (and hence, at rest 

relative to the rigid central carrier) , and with origin in the 
vicinity of the idealized attachment point (or points) 


X RB. y RB. Z RB. ~ the ( RB j)''f rame has origin at the CM of m RB and is 

oriented as the RB. -frame when m RB has not undergone a 

rotation relative to the carrier. (The (RBp* frame is at 
rest relative to m-,^, ) 

Kd . 


□ 


RB. 


1 = the inertia matrix of m RB referred to xJ^ B y RB z RB (context should 

l l i l 

make clear whether the symbol m-,^ denotes the body whose mass is 

KJj j 

or the numerical value of that mass) 


m 


RB. 


^PRB ’ ^YRB ’ ^rRB ~ Euler an gl es defining the orientation of the (RBp* 
1 1 1 frame relative to the RB.- frame 

l 


T rb = T (B+RIb) — a rotation matrix defining the transformation from the B 
1 resolution to the RB. resolution 

T' = T(RBi -v (RB-)’) - the rotation matrix defining the transformation from 

XVl3 • 1 

1 the RBi resolution to the (RBj)' resolution 

T k B . = t 1 '»yRB.'(3) [ ^PRB. 1 <2) for 8 2 ' 3 > 1 SeqUenCe thr0U S h 

^PRB- ’ *yRB.’ ^rRB. in P assin S from the RB i orientation 
to the (RBj)' orientation 

^RB. = T k Bi T RB. = T (§ - <RB i ),:i 

position, vector referred to the RB. frame, of the CM of m RR (the 

origin of the (RB^)' frame) when m RR has not experienced either 

i 

translation or rotation relative to the central carrier 

Are = displacement, relative to the carrier, of the CM of m RB . (This vector 
1 has the RB. resolution) 


r<°) = 

RB. 
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Jr B = angular velocity of the (RB.)' frame relative to the RB. frame (this 
i 11 

* — - A. T j. 1_ / T'* n x t * . • 
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i ill 


^RB v H “RB x v for any vector v having the (RB.)' resolution 
i i 1 


CMRB. 

R g 1 = position vector of the CM of m RB referred to the S -frame 
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^RB. = P° sition vector of a generic point of m RB referred to x^ B y' B z' B 
i i l l i 
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A rb = surface area of m RB 


= T^ b ^„ b = B resolution of w RB (also the B resolution) 

13 i i 1 i 


RB 
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*t sT ■!+ 

J RB. " 1 RB. W RB. 
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gRB. = torc l ue on m RB. alone due to Earth's gravity field 


MgR^-) _ torque on m RB alone due to Sun's gravity field 
-> (MOON) 

MgRB. = tor< ^ ue on m RB. 810116 due to Moon ' s gravity field 

(ATTACHMENT) _ gum ^he damping moment and restoring moment resist- 
1 ing, through the attachment, any rotation of m RB 

relative to the carrier . 1 


The equation governing the translational motion (relative to the central carrier) of 


m RB * s f° und to 

l 


— ^ 

C xRBi 

0 

0 ‘ 


K xRBi 

0 

0 “ 

m RBi A RBi + 

0 

C yRBi 

0 

A RBi + 

0 

K yRBi 

0 


. 0 

0 

C zRBi^ 


. o 

0 

K zRBi J 


1 RBi 


+ m RBi T RBi I 2 T RBi A RBi + ^ + ) [r RBi + T RBi ( *RBi + A RBi )] 


+ T 


R + 


'M s 


5 (MOON) ,2 (MOON) £ CMRBi, £ (SUN) 

K o Hm ( K o " K c 1 y s K s 


s |g (MOON) 3 
s 


W M Vil s -s 

g (MOON) j* CMRBi 3 
s s 


R. 


(SUN) ,3 


s 


,2 (SUN) 2 CMRBi, 

U (rv ~ K ) 

K S v S S 


^(SUN) £ CMRBi 3 
R - R g 


>- F - T 

gRBi x RBi 


f (dF AERO*(B) 


(A RBi ) 


- T 


RBi 


/ 


(dF SOLAR' ) (B) 96 0 ’ n rb 


(4- 94) 


(A RBi } 


80 



Clearly evident in equation (4-94) is the assumption of viscous damping and linear 

— >• 

restoring forces at the attachment. The symbol F denotes the RBi resolution 

git r>l 

of the Earth’s gravitational force on m R g. and is obtainable by an equation analogous 
to (A- 30). The procedure for finding should be clear after reading the 

paragraph containing equation (A- 30). In fact, that portion of the paragraph starting 

— >■ 

with its second sentence can be made applicable to in the following way: 

replace i by RBi when i is used as a subscript or superscript; replace the expression 

"flexible appendage i" with the symbol m RB ^ and delete the expression "in its deformed 

state;" delete also the remark within parentheses appearing as the last part of the 

final sentence of that paragraph; and ignore the references to Sections 2 and 3. In 

passing, one may observe that if m R g. is regarded as a point mass (having a null 

inertia matrix and zero surface area), if T^. = 0, if the two components A and 

^ Kol yKJal 

A zRBi of the vector A RBi 5 [A xRBi’ A yRBi’ A zRBi ] are ^trained to be zero, 

and if both members of equation (4-94) are premultiplied by ii^. = [1, 0, 

KJjI KJdI 

the result is an equation completely equivalent to (4-74), in content, that is, which 
is as it should be. 


The contribution of the relative motion of the m DD ., i = 1, 


Nr R , to the 


"system" translational motion is realized by adding to the right member of equation 
(4-53) the expression (4-95) below, and the contribution to the "system" rotational 
motion is expressed by the addition of expression (4-96) to the right member of 
equation (4-73) 
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If equation (4-73.14) is designated the "system" rotational equation, then the expres- 
sion (4-96) with r CM deleted should be added to the right member of equation (4-73.14). 

The second sum in expression (4-96) is equivalent to the expression (4-96A), 
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RB 


-£ < ? Li D RBi * 


+ (Wtj + q 


) x n RBi T 

) D T3 1 - II 1 D D.' 


RBi RBi T VW B T RBi' A RBi 1 - 1 RBi RBi 


w' •} , (4-96A) 


i=l 


which expression one could rightly expect after an inspection of equation (4-73). 
The rotational motion of relative to the carrier is determined by 

0 1 ( ^RBi + ^RBi W B + fi RBi ^RBi W B ) + (w RBi + ^RB. a> B ) 


x D 1 (a) RBi + ^RBi U B ) " ^gRBi 
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RBi 


(A RBi ) 


X ^RBi ^ dR SOLAR^B " f r RBi X ^RBi (d *AER(pB 

(A RBp 

_ (SUN) _ (MOON) _ jjj (ATTACHMENT) _ -* 

gRBi gRBi 1V1 RBi ’ 


i = 1, ... N 


RB 


(4-97) 


If the dimensions of m RB ^ are not sufficiently large to warrant retention of the gravity 

terms, solar terms, and aerodynamic terms in equation (4-97), as is likely to be the 

case, then they may be ignored along with terms arising from certain other sources of 

excitation mentioned near the beginning of this section (see the paragraph following 

equation (4-1G)). Even if deemed important, the gravity terms due to Sun and Moon 

should be adequately approximated via equations similar to (4-55) and (4-56). Equa- 

— ► f g UN } 

tions (4-55) and (4-56) are made directly applicable to the computation of M g RB j 
and by mer ely replacing the subscript B by RBi, the symbol Q by Q RB *, 

and the superscript CM in the definitions of Rg^ and Bjyjqqn by CMRBi. Just as 
FgRBi may be approximated by an equation similar to (A- 30), as explained above, 
so can be approximated by an equation similar to (A-31). 
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Not only would the presence of the m^g.., i = 1, ... Ngg , require the addi- 
tions to equations (4-53) and (4-73) as described above, but would result in slight 
change in the expressions for r^ , r^j , and [J as well. To the expression within 
braces in equation (D-l), one should add 
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RB 


i=l 


m RBi T RBi A RBi 


(4-98) 


while to that within braces in equation (D-2), one should add 
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RB 


m (r + ) 

RBi vr RBi A RBi ^RBr 


i=l 


Obviously, the necessary modifications to (D-3) and (D-4) are made by adding the 
first time derivative of expression (4-98) to the expression within braces in equation 
(D-3) and the second time derivative of expression (4-98) to the expression within 
braces in equation (D-4). Equation (E-4), defining the system inertia matrix [] , 
would be modified by the addition to its right member of the expression (4-99). 
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The vector q DD ., not defined in Appendix E, denotes the instantaneous position of 
the CM of m„ R . relative to the B -frame and is given by 


*RBi 


= r 


RBi 


CM 


+ T 


(it 


(o) 


RBi v RBi 


a rbP 


The subsequent modification of equation (E-5) is effected by simply adding to its 
right member the first time derivative of expression (4-99). 

No attempt will be made in this paper to assess the effect of liquid sloshing on 
vehicle motion since it is hardly possible that the author could add anything to the 
extensive literature on that subject. However, with regard to slosh, Reference 21 
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should be cited as one providing a collection of formulae pertinent to a variety of 
tank configurations plus an exhaustive list of technical papers on sloshing in both 
high-g and low-g environments. Due consideration cannot be given to slosh until 
tank geometry and mechanical analog of the liquid have been decided. Only then 
can one write the sloshing equations and add the appropriate terms to equations 
(4-53), (4-73), and (D-l). Introduction of the mechanical analog requires that the 
contribution of the fluid to the system inertia matrix [] (see equation (E-4)) is 
understood to be the "effective" inertia matrix of the fluid. The structure of some 
slosh models is such that certain of the equations (4-74) or (4-76) may be designated 
as slosh mass equations. 


V. REMARKS ON SIMULATION 

To digitally simulate vehicle motion, such a simulation being nothing more than 
the numerical solution of the system of differential equations descriptive of that 
motion, one must first recast the underlying equations in a form suitable to pro- 
gramming, that is, a form amenable to the direct application of the selected numerical 
integration scheme. Use of the popular fourth order Runge-Kutta technique, for 
example, would require that the system of equations be put in its equivalent first 
order form. 

Also necessary to the construction of a simulation program, as clearly indicated 
by the equations in Section IV, are many subsidiary relations, not the least important 
of which are the control equations. Linked to the control equations are expressions 
for approximating measurements made by onboard sensors. Equations for determining 
the output of the filter networks designed to filter those measurements should be 
implemented. It is hardly necessary to point out that the equations in the previous 
sections and the appendices are indicative of program input with the exception of 
such data as the time interval of interest, integration step size, etc. 

The reader who is concerned about the existence and uniqueness of solutions 
is hereby advised to consult the literature on ordinary differential equations. 

A proposed sequel to this paper is one wherein the equations of this paper are 
applied to one or more specific station configurations. 
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APPENDIX A 


ON GRAVITY FORCE AND GRAVITY TORQUE 

The manipulations leading to the expressions below for the force and moment 
exerted by the Earth's external gravitational field upon a body of arbitrary shape are 
based upon the gravitational potential function of Reference 6 which includes spherical 
harmonics through fourth order and reads, in the notation of this paper, as follows: 


U g 


( 4 n 
TT I 1 - E E 

^ n=2 m=0 


nm 


n 


= 1 P m (sin 6) cos m (A - 


n 


A ) 

nm 


(A- 1) 


wherein denotes the product of the mass of the Earth and the universal gravita- 
tional constant ; R , A , and 6 are the spherical coordinates of the field point relative 
to (see Sections 2 and 3); A nm is the longitude (positive east of the prime 

meridian) of the principal meridian of symmetry for the nm harmonic; the J nm » n = 2, 
3, 4, m = 0, ... n, are dimensionless coefficients peculiar to the planet Earth with 
J 21 - 0 by virtue of the assumption that the Zg axis (and also the Z g axis) is a 
principal axis of inertia for the Earth; a is the mean equatorial radius; and P™ (sin 6) 
is the associated Legendre function of the first kind of degree n and order m defined 
by 


,m 


P™ ( sin S ) = cos m 6 


P (sin 6 ) 
n 


d (sin 6) 


m 


the function P n (sin6) being the Legendre polynomial of degree n in the argument 
sin 6 . 


As seen in equation (A-l) the dimension of U is that of work per unit mass or 
potential energy per unit mass, the equivalent of the product of unit of force and 
unit of length divided by the unit of mass. From dimensional considerations alone it 
is evident that a differentiation of U with respect to distance in a given direction 

o 

(the directional derivative, that is) gives the gravitational force per unit mass, more 
frequently called the gravitational acceleration, in that direction. It can also be 
argued (Reference 7) that a differentiation with respect to arc length along the curvi- 
linear coordinate curves whose unit tangent vectors are the u^, u ^ and u ^ defined 
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in Section 3 gives the components in the directions of those unit vectors of A , the 
symbol A here denoting the acceleration due to the Earth's gravity at the field point 

o 

with spherical coordinates (R , A, 6) referred to the E-frame. Thus, 


A g E A gR 5 R + A gX + A gS U « 



(A-2) 


The results of performing the indicated differentiations in equation (A-2) are given in 
Reference 6 and are repeated here (in the notation of this paper) for ready reference. 
(The author of this paper has verified, through tedius scratchwork not shown, that 
the results appearing in Reference 6 are, apart from some obvious typographical 
errors, correct. Minor typographical errors are also present in the expression for 
the potential function in the reference cited.) 
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(A-3) 


90 



-i 9 U y 

1 _g _ 


gA R cos 6 9 A 
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+ | (g j J 31 (5 sin 2 6-1) sin (A - A 31 ) 


+ 30 ( J 32 cos 6 sin 6 sin 2 (A - A 32 ) 
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+ 15 | R J J 32 (3 sin* 6-1) cos 6 cos 2 (A - A 32 ) 


+ 45 ( R j J 33 cos 2 6 sin 6 cos 3 (A - A 33 ) 


I (!) J 4 Q (7 sin 2 6-3) sin 6 cos 6 


(A-4) 


2 - ^22^ 


(A- 5) 
(Continued) 
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(A- 5) 
(Concluded) 


+ 


+ 


+ 


+ 



41 (28 sin 4 6 -27 sin 2 6+3) cos (A - 

2 

42 (7 sin 6-4) cos 6 sin 6 cos 2 (A - A 

9 2 

J 43 (4 sin 6-1) cos 6 cos 3 (A - X 4 g) 
J 44 cos 3 6 sin 6 cos 4 (A - A 44 ) j , 


42 ' 


The force and torque* (attributed to the Earth's gravity field) experienced by 
a space vehicle** will here be denoted by F g and M^g, respectively, the second 
subscript B being indicative of the B -resolution. 


Definitions of symbols incident to the development of expressions for F^g and 
Mgg are in order now (some being repeated elsewhere). 

-* CM 

R = position referred to the S- frame (see Section 2) of the CM of the entire 

b 

vehicle system. 

r = xi + yj+zk = position referred to the B -frame (see Section 2) of a 

generic point of the vehicle. 

-» -* CM T -» 

R = R +T r = position of a generic point of the vehicle relative to the 
s s 

S -frame (T = T(S-*B) defined in Section 3). 

R Q = T R s CM = B -resolution of Rg CM , R Q = |R Q | = |R g CM | 

R=TR g = R o + r = B-resolution of R g , R = |R | = |R q + r| 

( u r) b = R/R = B-resolution of u R (see Section 3) 

Uro = R q /R 0 = B-resolution of R g CM /lR g CM | 


* The moment reference point is the vehicle center of mass (CM). 

** The body need not be a space vehicle but may be any body of arbitrary shape 
and finite dimensions occupying a position in the Earth's external field. 
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n = T K = B- resolution of K 
s s 




I 

i 

XX 

xy 

xz 


i 

I 

yx 

yy 

yz 

zx 

*zy 

i 

zz 


Inertia matrix of entire system referred to xyz (the 
B -frame with origin at the system CM). 


*xx = 

/< y 2 < 

- z 2 ) dm , 

‘yy = / (x2 + z2) dra - 

= (x 2 + y 2 ) dm 


(m) 


(m) 

(m) 

1— 1 

vn 

II 

- f 

dm , 5 ± 

n , = x, y, z 



(m) 




/<• 

. . ) dm = 

//A-* 

a (x, y, z) dxdydz , (a(x, y, 

z) = mass density) 

On) 


(m) 




(as used immediately above the subscript (m) indicates that the integration is to be 
extended over the volume occupied by the entire vehicle system. Elsewhere the 
symbol m will denote the numerical value of the mass of the entire system.) 


Tr(fl) = trace of D = 1 + I + I 

u L-l xx yy zz 


In approximating M^g the contribution of the longitude dependent terms (the 
tesseral and sectorial harmonics) was completely ignored while in approximating F^g 
their effect was accounted for "in part," that is, their variation over the vehicle was 
ignored, it being assumed that their values at any point of the vehicle differ 
negligibly from those at the vehicle CM. 


The contribution to F _ of the dominant term in the expression for A _ , here 
, ( n n \ 
denoted by ' , is given by 


F (o,o) _ 

gR 


■ f 


dF 


(m) 


gR 


(o,o) _ 


*E y^gr ? , 3 

(m) > R o + r l 


(R/-. + ?) 


dm 


m y. 


% 


u 


Ro 


3 y i 


R 


□* 


Ro 


3 y, 
2 


Tr® - • & DSroISro.CA-S) 

o L 
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while its contribution to M 


gB 


is 


M 


gR 


(o,o) 


/ 

(m) 


r x dF 


gR 


(o,o) 


3 y ■ 


% 


R 


o 



(A- 7) 


In arriving at the approximations (A-6) and (A- 7), use was made of the approximate 
expansion 


R 


K = IR 


+ r 


-K 


R 


-K 


- K R 


-K-l 


K R 


-K-2 


u 


Ro 


. K (K+2) 
2 




-> -+T 


r r 



(A-8) 


and the relations 


/ 


r dm = 0 


(by definition of the mass center) 


(A- 9) 


y r T r dm =| 7r([]) 

(m) 

f * ? T dm = - [|+ i 7^(0) I (3x3 ) 
(m) 


*(3x3) 



0 

1. 

0 


0 

0 

1 . 


(A- 10) 


(A- 11) 


ilk 

in addition to further simplifications made by discarding all terms in x y J z when 

5 

i + j + k > 3 (such terms have denominators of the order of R q and greater). 


It should be remarked that the axes xyz (the B-axes) have not been restricted 
to be principal axes of inertia and further that the inertia matrix [] must be inter- 
preted as that of the instantaneous deformed configuration referred to the xyz axes 
which have origin at the instantaneous position of the CM. 


In finding the contributions of the zonal harmonics to F _ and M^g (these 
contributions being identified herein as M gR^ 0 ^’ F g6^ k, °^’ ^g6^ °^ 

k = 2, 3, 4) the relations (A-8), (A-9), (A-10), and (A-ll) again find repeated 
application as do the following: 
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C -'5- 



^ Ip ^ 

sin 6 o = ^Ro n ’ ( 6 o = S eocentric latitude of vehicle CM) , (A- 12) 


sin 6 n T (u R ) B R | 

f ^ r n v \ 

I s ” + R ) 

(A- 13) 


^ o / 


cos 6 (vf^ ) B = n - sin 6 

^R^B 

(A-14) 


the relation (A- 14) being invoked only in the development of expressions for the 
^^,5 ^ ^ anc l k = 2, 3, 4. Notice that the use of (A-14) circumvents the 

need for an approximate expansion for cos 6 . 


No exposition of the detailed manipulations leading to the approximations below 
will be given, it being presumed that the reader can, with the time and inclination, 
supply all the manipulations omitted. 


3 a 4 J 9n r (3 sin 2 6-1) (3 + r) 


? (2,o) _ r ^ (2,o) _ y E a °20 r ' ' o 

gR */ gR 2 J 


dm 


(m) 


(m) 


% J2 ° { (3 sin 2 V l) + (7 sin 2 s -2) T t( []) 

R 

o o 


+ 2R 7 <35 ' 189 d " 2 6o> “R° 


Ro 


R 


n T [}n 


42 sin 6 


R 


o -+t n- i -* 

- n U u Ro } u Ro 


3 y E a J 20 
2 R 6 


{(21 sin 2 5) IQ - i 7J(j], I (3x3)] iI Ro 


- 6 sin 6 o I D - 2 A(0) I (3x3) ] S} 


(A-15) 
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M 


(2,o) _ 


■ f 


r x dF 


(2,o) _ _ 3 y E a J 20 ( - f (3 sin 2 6-l)r 


gR ‘ ~ gR 

(m) 


{ S o x f 


dm 


(m) 


R' 


% - 


3 y E a J 20 
2 R 5 


{(21 sin 2 6 q -5) u Rq x (Qu R o ) - 6 sin 6 q u Rq x c[]H) 


+ 3 sin 6 q 7r([|) (u Rq x n) } 


(A- 16) 


gR 


(3,o) 


- /'av (3,0) Q \j T 

= J dF_ t( v = 2 P r a J 


(5 sin 6-3) sin 6 (R + r) 


(m) 


gR 


3 - / ° 


E 30 


(m) 


R 


dm 


^ 2 y E a3 J 3o { rr (5 sin3 6 o ' 3 sin V 

R 


+ — (15 sin 3 6 q - 8 sin 6 Q ) Zr([]) 
R 

o 


2R 


(21 sin S o - 55 sin 3 6 q ) u Rq []u 


Ro 


15 sin 6 


R 


2 n T 0" + A CW5 sin 2 « 0 -21) H T []u Ro ) u Ro 
R 


2 y E a J 30 


R 


T~^ { 3 sin a o <15 sin 2 « 0 -7)[[|- \ AcQ) I (3x3) ] H Ro 


+ 3<1- 5sm 2 6 o )[Q^ 1 7i<Q> l (3x3) ]H) 


(A- 17) 
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M 


gR 


< 3 -°> = /?* d? gR <3 -°> = - 2 u E a 3 J, # { R 0 x - 3 sin 6)? dm J 

(m) (m) R 


% 


-2 p £ a J 3Q 
R 6 


< 3 sin a o (15 sin 2 « Q - 7) u R(J x ([]3 Ro ) 


3 (1-5 sin 2 6 q ) u Rq x (Qrt) - § (1-5 sin 2 6 Q ) 7"r(|]) (S Ro x n) } , 

(A- 18) 


p (4,o) /" 

gR J 


dF 


(m) 


gR 


(4,o) _ 5 y E a J 40 


y' ~Y (35 sin 

r__x R 


'- 4 5 - 30 sin 2 6+3) (R Q +r) dm 


(m) 


% 1^40 | . (35 sin 4 6o 

K 


30 sin 2 6 q +3) 


+ — ig (385 sin 4 6 q - 300 sin 2 6 + 27) 7r([|) 

R o 

- (5005 sin 4 6^ - 2970 sin 2 6 + 189) uj (1^ 

2 j^ ° o o Ro U Ro 


+ ( 1-7 sin 2 S Q ) n T 0S + j/ 0 (77 sin 2 i Q - 27) n T |]u Ro } 5 

o o 


Ro 


5 a J.q (A 9 n 

+ s < (385 sin 4 6 - 270 sin 6 +21) II u„ 

sr 8 \ ° o u r 0 


- 20 sin 6 Q (7 sin 2 6 Q -3)[]n - ^ (385 sin 4 6 q - 270 sin 2 6 Q +21) 7r([]) u 


Ro 


+ 10 sin 6 Q (7 sin 2 6 Q -3) 7r([j) n j 


(A- 19) 
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(4,0) _ 


M = / ? x df gR <4 ' 0> = - SVB » J - P { 5 o x fh < 35 


gR 


(m) 


. 4 

sin <5 


(m) 


R 


- 30 sin 2 6 +3) r dm j 


0 j 


~ 5 1J E a J 4o 
7 

8 R ' 


{(385 sin 4 <5 q - 270 sin 2 6 q +21) u Rq x (Qu R o ) 


- 20 (7 sin 3 6 o - 3 sin 6 Q ) u Rq x ([]n) 

+ 10 (7 sin 3 6 q - 3 sin <5 Q ) 7r([]) (u Rq x ri) } 


(A- 20) 


g<5 


(2,o) _ 


= y'd? 


(m) 




(2,o) _ 


= -3 a J 




E “ "20 y „4 

(m) R 


sin 6 cos 6 (u.) D dm 

0 D 


- 3 P R a J 2Q . 


{ y^- 


dm n + 3 p n a J 


E “ 20 


(m) 


/ 

(m) 


sin 2 6 


R 


5 v o 


(R + r) dm 


% 


3 a J 20 
R 4 


{ si " 2 S o[ 


m + 


[ 7r(P) 

2R 2 2R 2 ~ R ° 

o o 


63 jT 


0 




14 sin <5 


R 


o 


3 a J 20 
R 4 


I sin 6 m + 
{ o 


TrQ) 


V 2R 


35 -*t n ^ ~1 , 5 ->t n -*• \ -* 

^"2 U Ro U u Ro + ^ n U u Ro }n 
o Jo 


+ 3 “e I , sin 2 So Qe!ro . , sin 4 


R 




■]i 


, (A-21) 
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( 2 , 0 ) _ r t v .dr (2,o) _ 


V <2 ' 0, = 3 »e> 2j 2.-{/tt ? dm } 

(m) R 


(m) 


3 y E a J 20 R 


„*{/ 


(m) 


sin 2 6 
R 5 


r dm 


} 


% 


3 y E a J 20 
R 5 


| n x (Qn) + sin 6 Q [7 sin 6 q u Rq x (Qu R o ) 
- 5 n x (|]u Ro ) - 2 u Rq x (|]n) 

+ | Tv<\\) (it x u Ro )] } , 


(A-22) 


3 (3,o) _ 

F g6 


/ dF (3 ’ o) - 


~ 3 y E a J 30 


(m) 




y* -Tj (5 sin 2 6-1) cos 6 (u^ ) R dm 
R 


(m) 


■{ 


~ 3 y E a J 30 


y -^g (5 sin 2 6-1) dm|n 
r \ R 


(m) 


3 )Jp a Jon 1 O 

+ g- - J (5 s in 6 6 - sin 6) (R q + r) dm 

(m) R 


% 


~ 3 y E a J 30 


2R 


g {m (5 sin 6 Q -1) (n - sin 6 q u R q ) } 


o 


3 yg a J 


2R 


{[5(1-7 sin 2 S Q ) TrcO) * f (» s* 2 «„-» 3^ []u Ro 


5 n T []n - 70 sin 6 q rt T []u Ro ]n + [3 sin 6 q (15 sin 2 6 q - 1) 7r([] ) 


+ | sin 6 q (7-55 sin 2 6 q ) u Rq Qu Ro - 15 sin 6 q it T \]n 

(A-23) 

+ (135 sin 2 6 0 -7)n T []u Ro ] u Rq + sin 6 q (45 sin 2 6 Q - 7) D^ Ro (C ° nt ’ ) 
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+ (1-15 sin 2 6 o ) □ n + \ [sin 6 q (7-45 sin 2 6 q ) u Ro + (15 sin 2 6 Q -l)n ] /r(Q) j , 


(A- 23) 
(Cone.) 


M 


g<$ 


(3,o) _ 


y?. *<'■«>. s* { 


"E 8 J 30 y J_ (5 sjn 2 j.jjJ dm J 


(m) 


On) 


- R 


o x { 


3 y E a J 30 


y -^r (5 sin 3 6 - sin 6 )r dm j 
, R 


On) 


15 y E a J 30 ( . 2 , 

% g { (7 sin 6 


2R 


3 y E a J 30 
2R 6 


{ < 7 sin2 6 o" 1) [" x (D 3 RO ) ' \ 7 r( D ) X ^Ro ) ] 

- 2 sin <$ o n x (Qn) j 
{ sin { o <45 s ” 2 6 o" 7) a Ro x <O a Ro ) 

o> [ S Ro x <0 3 > - \ 7^<D) < a Ro x B >] } ’ 


+ (1-15 sin 2 6 


(A-24) 


(4,o) 


dF 


g« - y “‘gs 

(m) 


(4,o) _ 


5 y E a J 40 


y ' (7 sin 2 6-3) sin 6 cos 6 (u^ ) fi dm 

■ v R 


(m) 


■{ 


-5 y g a J 4Q 


y (7 sin 2 6-3) sin 6 dm | n 
On) R 


5 y E a J 40 


y -y (7 sin 2 6-3) sin 2 6 (R q + r) dm 

\ R 


(m) 


% 


o 21 sin 6 o - 7 - i— i 

| m (7 sin 3 6 q -3 sin 6 Q ) + 2 ( 3 s ^ n ' r (|J) 


2R 


R„ 


(A- 25) 
(Cont. ) 
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63 sin <5 o T n 21 sin 6 


“ <3-11 sin 2 5 0 ) u Ro 0“ 


R° 2 

o 


o ->t n -► 

- n U n 


(9sin 2 6 0 -l) n T []3 Ro } n 


5 P E a J 40 


| m (7 sin 4 <S Q - 3 sin 2 6 q ) 


+ — ^ (154 sin 4 6 q -39 sin 2 6 q -3) 7r([]) 
2R 

o 


• 2 , 
sin 6 


- (297-1001 sin 2 6 q ) u Rq \}Z Rq + (1-14 sin 2 6 q ) n T \]n 


2R 2 — ““ o' “Ro LTRo T ^"2 

O o 

2 sin <5 

+ —TT 2 - < 154 sin V 27 > S 0 a Ro} 3 Ro 


5 P E a J 40 ( .2 


{ sin 2 « D (77 sin 2 V 27 > [Qu Rq - 4 ^(Q) u R q ] 


+ 2 sin <5 o (3-14 sin 2 < 5 q ) [[]n - ^ 7r([]) n] J 


(A- 25) 
(Cone. ) 
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M 




<4 ’ 0) = 


5 y E a J 40 


(m) 


-^g (7 sin 2 6-3) sin 6 r dm 
(m) R 


\ x { 


5 u E a J 4Q 


y -y (7 sin 2 6-3) sin 2 6 ? dm j 

' \ R 


(m) 


% 


5 Pg 8 J 


2R 


^ | 21 sin 6 q (3 sin <$ 0 - l) n x [[]u Rq ~ 2 ) u Ro^ 


-2 sin 6 q (3-14 sin 2 6 Q ) u Rq x [Qn - ^ 7r(|J) n] 

+ 3 (1-7 sin 2 6 q ) n x ([] r?) - sin 2 6 Q (77 sin 2 6 Q -27) u Rq x (Qu R o ) } . 

(A- 26) 

As remarked in a previous paragraph , the variation of the longitude dependent 

— ► 

terms (in the expressions for the components of A ) over the vehicle will be ignored, 

o 

it being assumed that their values at any point of the vehicle are equal to their 

values at the vehicle CM. An immediate consequence of this assumption is that the 

approximation to M^ R is devoid of any contribution from the tesseral harmonics whose 

contribution to F D must then be 

gB 


? gB = mT I“ol?3) < S g >(o) 


(A- 27) 


In equation (A- 27), m denotes the mass of the entire vehicle system; T = T(S^-B) is 
the rotation matrix defined in Section 3; the angles a Q and 6 q are, respectively, the 
right ascension and declination (geocentric latitude) of the vehicle CM, these being 
given by (among other expressions) 


, -1 , v CM /v CM. „ „ 

a Q = tan (Y g /X g ) , 0 < a < 2tt , 


. . -1 CM /13 CM. /0 . 

6 = sin (Z /R ) , - tt/ 2 < 6^ < ir/2 , 

U 00 o 
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where X Y ^ and Z are the components of R ^ and R CM = 

s s s s s 

A* denotes the "longitude dependent part" of with u^, and 

defined by 


R CM |; and 
s 


components 


* 

gR 

= A „ 

gR 

minus 

the 

terms 

in J 20 ’ 

J 30’ 

J 4 0 and minus the term 

* 

g* 

= A , 
g* 







* 

= A , 
g<5 

minus 

the 

terms 

in J 2Q , 

J 30’ 

and J 4Q , 


E 


“** i CM i 

the subscript (o) on A^ indicating evaluation at R = R q = |R g |, <5 = S Q , X = X Q , 

where X , the east longitude of the vehicle CM, is given by 


A o = a o ' °V " “il (t ' V ’ 0 < < 2 tt 

From the foregoing, the approximations to F 
this paper, read as follows 


gB 


and M 


gB’ 


in the notation of 


F % F * 
gB * gB 


♦ L 

K=0, 2, 3, 4 


$ (K,o) 

V 


+ 


z 


8 (K,o) 

F g« 


(A- 28) 


M 


gB 


% 


M 


K=0, 2, 3, 4 


gR 


(K,o) 


K=2,3,4 


M 


g6 


(K,o) 


(A- 29) 


The reader should be aware that the dependence of F^g and M^g upon vehicle attitude 
enters through the dependence of the unit vectors u Rq and n upon the Euler angles 
specifying the orientation of the B-frame relative to the S-frame. One should notice 
also that it has not been necessary to introduce an additional reference frame (an 
"orbital" frame whose attitude relative to the S-frame also changes); and therefore, 
no additional attitude reference angles, thereby circumventing the need to transform 
the inertia matrix. The more prominent writers on the subject of gravity torque have 
used an orbital reference system (such as that alluded to) in defining "attitude devia- 
tion" angles in terms of which they expressed the potential function (via the trans- 
formed inertia matrix) and subsequently differentiated the potential function with 
respect to these angles to get the components of gravity torque. 
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Though the expressions for F _ and M D were derived with the entire vehicle 

grs gij 

system in mind, they are applicable to any component (as if the component were up 
there alone) by a proper interpretation of the symbols, that is, by an appropriate 
alteration and/or change in the meaning of the symbols. In particular, they can be 

+ Vi 

made to apply to the 1 flexible appendage of the space station model described in 

Section 1, and use will be made of them in writing the contribution of gravity to the 

1 h 

equation corresponding to the coordinate 0. (associated with the i- flexible appendage 
and denoting the rotation of that appendage, as a whole, relative to the rigid central 
carrier) . 

The symbol F ., denoting the i-resolution (see Section 2 for the definition of 
the x.y.z. frame) of the Earth's gravitational force on appendage i alone , presents** 
itself in the equation corresponding to the coordinate 6.. On introducing the symbols 


rm = mass of flexible appendage i 


(A gi } (o) = K evaluated at R = R Q ., X = A^, 6 = S Qi 


rt CMi _ rv CMi v CMi „ CMi.T ... e , , ,, „ , „ ,, 

H s - [X g , Yg , Z ] = position referred to the S-frame of the CM 

of flexible appendage i (in its deformed 

state) 


- 1 CMi CMi 

a oi = tan ( Y c /X ) = right ascension of the CM of the flexible appendage i 


(0 < a . < 2ir) 


0 i “ s ^ n '(Z^i/lR CMi|) _ g j n i(n T u ) = declination or geocentric latitude 

of the CM of appendage i 


(- tt/2 < 6 Q . < tt/2) 


- “ 0 i " “p ‘ aj ^ (t - t Q ) = east longitude of the CM of appendage i 


(0 < X o . < 2 tt) 


R q - = T^ T R g CIVR = i - resolution of R (see Section 3) 


R oi = R oi = R s 


CMi 


u Roi *W R oi 


= T\ T K g = i-resolution of K g (see Sections 2 and 3), 
104 **As one factor of a scalar triple product. 



it should be obvious, without recourse to the method of deriving (A- 28), that 


F . & 
gi * 


T. T 
1 


[ot oi ] ( 3) 


l- S oi>?2) 


I m i 


(A*. ) 


gi (o)j 


+ 


X 

K=0 f 2, 3 , 4 


S (K,o) 
F gRi 


♦ E 


a (K,o) 
F gSi 


(A- 30) 


the F^ R . V ’ and F gfii ’ * n e Q uat i° n (A- 30) being found by simply replacing m, 

R o , u R q , n, sin <5 q ,Q , and Tr([] ) in the expressions for Fg R ^ k ’ 0 ^ an< ^ F g6^ °^ 

by nn, R q ., u R q ., rr, sin <$ Q ., Q *, and Tr([]*) 5 respectively. The symbol Q 1 is to 

be interpreted as the inertia matrix of flexible appendage i (in its deformed state) 
referred to the axes x!y!z! defined in Section 2 (recall that the axes x!y!z!, though 
oriented as x^y.z^, have origin at the instantaneous CM of flexible appendage i) . 


Also, appearing (as one factor of the scalar product of two vectors) in the 

— ► 

equation corresponding to the coordinate 6. is the symbol M . representing the 

O 

i- resolution of the gravity torque on flexible appendage i alone about its instantaneous 
CM. By a direct application of equation (A- 29), one has immediately the approxima- 
tion 


M . a» Y M 
gl ^ LmJ gRi 


(K,o) 


+ 


K=0,2,3,4 


X M 

K=2, 3, 4 


(K,o) 


g6i 


(A-31) 


the expressions for the Mg R / k, °^ and b e i n g obtained from those for 

Mg R ^ k ’°^ and by merely replacing R q , u R q , n, sin 6 Q , [] and Tr(Q) by 

R oi> u Ro ., rn , sin and Tr([] 1 ), respectively. 

— ► — 

In the equation corresponding to the coordinate 0. the terms in F . and M . 

o o 

can be combined to give the expression 



{M . + (a. (o) + Y (l) 

gi l 


n (i) ) 


x F .} 
gi 


(A- 32) 
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which, after some reflection on the matter, one should recognize as the x^ component 
of the gravity torque on flexible appendage i alone about the origin of x.y.z., (see 
again the coordinate system definitions of Section 2). In the expression (A- 32), the 
vector &/ 0 ) is the position vector, referred to x.y^z., of the CM of appendage i when 

= 0, that is, when appendage i is undeformed; is the Nj x 1 column matrix 

of generalized bending displacement coordinates n^, ••• (functions of time to 

be determined); and the 3 x N. matrix ¥ ^ has for its j th column the column matrix 

defined by 


(i) 


= — /*?• 
m. J ] 


(i) 


dm 


j = 1. 


N i 


(A- 33) 


(mj) 


the 3x1 column matrix <p ^ being the j** 1 column of the undamped modal matrix 4> ^ 
associated with flexible appendage i, i = 1, ... NA. The problem of deciding the 

number N. , that is, an adequate number of mode shape functions V > . 1 to include in 

1 . (il J 

the truncated modal matrix $ , will not be addressed in this paper since suitable 

criteria for selecting the modes to be accounted for have been the subject of investi- 
gation by more learned men whose work will be found in the literature. It should be 
remarked (if not already clear to the reader) that it is herein supposed that 'P is 
a function of position referred to the i- frame so that one could write (in the usual 
functional notation) 


^. (l) = <P . (r.) = 
1 J i 


ffi. (x. ,y. ,z.) q 
]x v i ,J V l I 

tfi. (x.,y.,z.) 
jy l u’ i 

<P. (x.,y.,z.) 

]Z l •’l 1 


, j = 1, ... N. 
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APPENDIX B 


ON AERODYNAMIC FORCE AND AERODYNAMIC TORQUE 


At space station altitudes, the flow regime is, presumably, free molecular flow. 
Following the development in Reference 10 pertinent to such a regime, the aerodynamic 
force, in the notation of this paper, is given by 



rdF (x) i 

Qr AERO 

^AERO f 

dF (y) 

at AERO 

(a t > 

dF 

L AERO J 


(B-l) 


where the subscript A T on the integral sign indicates that the integration extends 
over the vehicle surface* and 


dF 


(?) _ 


AERO 


= (1/2) p a V R { t°t (E 1 y 1£, + e 3 y 35 } + (2 - a n^ e 2 Y 2^ 
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^2 (1 + erf e 2 S) + 

s /IT 


2 „2 
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2 ? 


exp (-e“ S") | + -^ (2 - a n ) (1 + erf e^S) 
2S 


°n y 2 ? 


T* r 

w 

T! 


£ 2 ^ 


g— (1 + erf e 2 S) 


+ \ exp (-e 2 S 2 ) 

s 


]} 


dA 


? = x, y, z 


(B - 2) 


The symbol p denotes local atmospheric density. V n is the magnitude of the relative 
a ->■ K 
velocity vector, V R , herein defined** by 


*The vehicle surface is here presumed a convex surface, thereby ruling out the 
effect of molecules reflected from other parts of the body upon the force on a 
differential element of surface area. 

**Note the tacit assumption that the velocity of all points of the surface is the same 
as that of the vehicle CM. 
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I 


-+ -4- fM 

V R = T (-R S UV1 ) , T = T (S->B) , (B-3) 

The e^, i = 1, 2, 3, are the direction cosines of the relative velocity vector in the 
local coordinate system, the local y-axis being directed as the inward normal to the 
surface and the local x and z axes tangent to the surface with arbitrary (but 
sensibly chosen) directions. If T RR = T (B+LOCAL) denotes the rotation matrix 
defining the transformation from the B resolution of a vector to its resolution along 
the local axes, then the e. , i = 1, 2, 3, are determined by 



The triples y^, ? - x > y> z > are > respectively, the direction cosines of 

H ► ->• — > 

the unit vectors l, j, and k in the local coordinate system, that is, 


' Y lx' 

— »> 

■ n 


Y 2x 

= t bl 1 = T BL 

0 

(B-5) 

, Y 3x 


. 0, 


Y iy 

-> 

■ o' 


Y 2y 

= t bl ] = T BL 

1 

(B-6) 

. Y 3y, 


. 0, 



tsi 

r— 1 


‘o' 


->• 


Y 2z 

= t bl k = T BL 

0 

. Y 3z . 
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The symbols T' and T! denote, respectively, the surface temperature and tempera- 

W 1 

ture of the incident molecules. The dimensionless quantity S, called the molecular 
speed ratio, is defined by 



the symbol R denoting the appropriate gas constant. Current best estimates of the 
tangential and normal reflection coefficients, cr^ and o n , respectively, restrict both to 
the range 0.8 to 1.0, completely diffuse reflection corresponding to = a R = 1. 

These dimensionless coefficients are also called momentum accommodation coefficients 
and sometimes momentum exchange coefficients. 

The symbol erf (£2^) denotes the error function with argument e^S, that is, 

e 2 S _ x 2 

erf (egS) = - p J e dx 

/ TT 

o 

In a simulation program, the parameters p and T. T should be available from 

a l 

atmosphere tables based upon an acceptable atmosphere model while the assumption of 
a constant T ’ would not be unlikely. 

From equations (B-l) through (B-7), it should be obvious that F^ERO has t * ie 
B-resolution as does the aerodynamic moment, M^ERO’ given by equation (B-8). 



AERO 


/ r x dF AERQ 

(a t ) 



(A t ) 


dF 

dF 

dF 


(x) 
AERO 

(y) 

AERO 

(z) 
AERO 


(B-8) 


wherein r denotes the position referred to the B- frame of the point of application of 
dF 

AERO* 
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Except for plane and cylindrical surfaces , the surface integrals in equations 
(B-l) and (B-8) defy closed form evaluation. In the case of conical surfaces and 
spherical segments, certain simplifications are possible, but even after such simplifica- 
tions are made, there remain integrals which yield only to numerical methods. 


As pointed out in Reference 10, a body in free molecular flow does not alter 
the flow, thereby permitting one to subdivide the body into a finite number of "simple" 
bodies (flat plates, cylinders, cones, etc.), determine the contribution of each part 
to a nd , and then combine the individual contributions by simple addi- 

tion. 


In the reference cited above (Reference 10), the development of a general 
expression for the differential of force coefficient corresponding to a prescribed 
direction is followed by the development of expressions for the normal force 
coefficient , axial force coefficient , and moment coefficient pertinent to a flat plate , 
a circular cylinder, a right circular cone frustum, and a spherical segment, the 
developments being valid for complete diffuse reflection. In the notation of this 
paper, the differential of the force coefficient corresponding to the spacecraft* 
5-axis direction (5 = x, y, z) is obtainable from equation (B-2) by dividing both of 
its members by \ p Q V ^ ^REF denoting a reference area), setting = a n = 1, 

and replacing T ' by T ’ (the symbol T ’ denoting the temperature of the reflected 
molecules), the result being 


dC { = { <E 1 Y U + e 2 y 25 + e 3 y 35 ) e 2 (1 + erf e 2 S) 

Kiir L 


+ — — exp (-e 9 S*) + (1 + erf e 9 S) 

S 2 J 2S 2 2 


2 c 2 J ^ Y 2§ 


I [ 


e 2 ^ 


li ___ 

2 / T! L S 


(1 + erf e 2 S) + exp (-e 2 S 2 )] } dA , (B-9) 

s 


-]) 


? = X, y, z 


*One should be aware that the spacecraft 5-axis direction does not, in general, 
coincide with the "local" 5 -axis direction. 
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In Reference 10, the development of equation (B-9) precedes that of (B-2) 
which follows from (B-9) via introduction of the reflection coefficients and o , the 
purpose of the reflection coefficients being to admit both specular and diffuse reflec- 
tion (complete specular reflection is realized when = a n = 0) . 


Once the dimensionless coefficient C 
found via the familiar equation 


5 


is known, the force component F 


U) 

AERO 


is 


F (5) 

r AERO 


= (1/2) C ? 


p a V R A REF 


5 = x, y, z 


(B- 10) 


If and C^, % ± n> denote, respectively, the force coefficients pertinent to 
the B-frame £-axis and rraxis directions, then the differential, dC^, of the moment 
coefficient pertinent to the B-frame ?-axis, the <;-axis being orthogonal to both the 
5 - axis and the n-axis, is expressible in terms of the differentials dC^ and dC^. 
Subsequent integration then gives C^, which determines the ^-component of 
in accordance with 


M 


(?) 

AERO 


= ( 1 / 2 ) 


p a V R C £ A REF L REF 


(B- 11) 


the symbol L REF denoting a reference length. 

An alternate expression (see References 9 and 12) for F^ EE q, som ewhat more 
tractable than (and deemed "inferior" to) that above is the following: 


F O/ 

AERO * 


p a v l f « 2 - °n - V 

(A x ) 


1 ) 2 1 

rr n 


+ ° t (£ v * V V dA * 


(B- 12) 


In equation (B-12), both the unit vectors £ n and i are presumed expressed 
on the B- vector basis, the vector Sl n being the unit outward normal to the differential 
element of surface area dA and 1 defined by 


£ v = V V R 


> 


the in equation (B-13) being the negative of that defined by (B-3). 


(B- 13) 
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APPENDIX C 


ON THE VECTORS F gOLAR AND ^ SOLAR 

References 11 and 13 (among others) provide the following expression for the 
force attributed to direct solar radiation 


P SOLAR P SOLAR C (1 + C rs } COS ® + (2/3) C rd^ n 

(a t > 


+ [(1 - C ) sin 0] "r} cos 0 dA 

rb 


(C-l) 


In equation (C-l), the unit vectors n and \ define, respectively, the outward 

normal and tangential directions to the differential element of surface area dA; the 

angle 0 is the angle between the incident ray and n; the symbol C denotes the 

rs 

coefficient of specular reflection (the fraction of incident radiation reflected specularly) 
and C r( j is the coefficient of diffuse reflection (the fraction of incident radiation 
reflected diffusely). The unit vectors n and rare here presumed to have the 
B- resolution so that F gQLAR *°° ^ as B -resolution. The value assigned to 
PgOLAR’ the solar radiation pressure, depends upon the vehicle's position. If the 
vehicle lies within the Earth's umbra p gQ L AR = 0; if the vehicle is within the penum- 
bra, 0 < P S0LAR < (MAXIMUM SOLAR RADIATION PRESSURE IN VICINITY OF 
EARTH); and if the vehicle lies in neither umbra nor penumbra, the value assigned 
to P SOLAR * n a s i mu l a ti° n program is likely to be its value at 1 AU (one astronomical 
unit) . 


Regarding equation (C-l), it should be remarked that it does not account for 
the effect of radiation reflected from one vehicle part upon another part or shading 
of one part by another. Furthermore, it should be obvious that the integration need 
not extend over the entire surface area but only over that part exposed to direct 
solar radiation. 


Except for very special configurations, not likely to be encountered in practice, 

a closed form expression for the surface integral in (C-l) is out of the question. 

— 

However simple the surface configuration may be, the evaluation of FgQLAR w * u 
probably always rest upon numerical methods for evaluating integrals. 
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With r denoting (as usual) position referred to the B-frame, having origin at 
the vehicle CM (see Section 2), one can define the solar radiation torque about the 
system CM by (C-2). 


M SOLAR P SOLAR 


y r x {-[(1 + C rg ) cos $ + (2/3) C rd ] n 


(A t ) 


+ [(1 - C ) sin 0] x} cos § dA 

fS 


(C-2) 


Assuming it possible to subdivide the surface area into a finite number of 
"sub-areas" on the i tJl of which is exerted the resultant solar radiation force FgQLAR 
at the point with position vector r . relative to the B-frame, one can write M g0LAR 
as the sum of vector products 


^SOLAR " H r cpi X P SOLAR ’ (C ' 3) 

i 

In general, one can go no farther than the integral expressions of equations (C-l) 
and (C-2). 
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APPENDIX D 
ON THE VECTOR r 

cm 

To find an expression for ? cm > one has only to invoke the definition of the 
center of mass to write 


/ r dm = I (r - r ) dm = 0 , 

J cm 

(m) (m) 

subdivide m as indicated by the subscripts on the integral signs in equation (4-16), 
and integrate after substituting from (4-25), (4-28), (4-31), (4-34), (4-35), (4-37), 
(4-39), (4-41), (4-42), and (4-43), attention being paid to (4-52.1) and (4-52.2) in 
integrating over rm. The result is expressible as 


~ ~ (o) 1 
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CM CM m 
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~ (o) ~ T -*• 

where r^' , T. 0 and A Eio are given by 
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that 


~ T ~ T T 

T io = < T i >e-=o = T i 
1 


A Eio ^ A eP 3 .=3 .=0 

yi pi 


The integral in equation (D-2) is herein supposed a known function of time so 

(o) 4 . 

cm 


r is presumed a known function of time. 


Assuming the time derivatives of m negligible leads to the approximations 

NA 


? * ? (o) + I 
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APPENDIX E 


ON THE MATRICES [] , Q 1 AND I 1 

Among the symbols necessary to a "working" expression* for the inertia matrix 
Q, there are some which have not yet been defined. These include the following 
(with the exception of whose definition will be repeated for ready reference): 
q , the position, relative to the B- frame, of the CM of the rigid central carrier; 
q Q , the position vector, relative to the B-frame, of the CM of the rigid central 
carrier; Ip*,, the inertia matrix of the rigid central carrier referred to axes (oriented 
as both the B and B frames) with origin at its CM; I , the inertia matrix of the fluid 
referred to the B-frame (further remarks regarding I will be made later); and the 
vectors q., q £ ., q R ., q p ., q Gi , q^., q 0G , Gig ’ and ^g’ denotin g> respectively, the 
instantaneous positions relative to the B-frame of the CM of flexible appendage i, the 
CM of swivel engine i, the CM of rotor i, the i^ d point mass, the CM of the gimbal 
of the i th SDOF CMG, the CM of the gyro element of the i th SDOF CMG, the CM of 
the outer gimbal of the 2 DOF CMG in question, the CM of the inner gimbal of the 
2 DOF CMG in question, and the CM of the gyro element of the 2 DOF CMG in 
question (the words "in question" being used to adhere to a previous agreement that 
an additional subscript i, or superscript i, would be suppressed on all symbols relat- 
ing to a 2 DOF CMG). For a specific vehicle configuration, q will be a known con- 
stant vector and 1^ a known matrix of constants, while I 1 , though variable, "can be 
determined." By simply stating that 1^ is "susceptible to being determined," the 
author has tacitly avoided the construction of a mechanical system whose motion 
duplicates the response of the fluid. Introduction of a mechanical analog would, of 
course, increase the number of system coordinates and require that more terms be 
added to both the translational and rotational equations. Knowing the rate at which 

mass is depleted through propellant consumption and assuming negligible sloshing, one 
f 

might define I adequately via tables with time as argument. 

In terms of symbols already defined, the position vectors (referred to the 
B-frame) of the previous paragraph must be given by 

*By "working" expression is meant one which can be put to practical use such as 
serving as a guide to a programmer in coding a subroutine (of a simulation program) 
whose function is the construction of fl • 
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% q o ' r CM 


q 4 = r. - ? CM + T. (£. (o) + ¥ (l) n (l) ) , i = 1, ... NA , 


q Ei r Ei r CM £ Ei A Ei ’ 1 ~ lj **• NSE 


q Ri r Ri ‘ r CM ’ i - 1. ... NR , 


q Pi r Pi " r CM + 5 Pi A Pi ’ i - 1. ••• NP , 
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i = 1, ... NSDOF 


q gi r gi r CM ’ i ••• NSDOF , 


q OG r OG " r CM 


q IG r IG ' r CM » 


q g r g " r CM • 


CORRESPONDING TO EACH 2DOF CMG 


In deriving an equation for the computation of the inertia matrix I 1 (defined in 
Section 4 in the paragraph containing equation (4-69)), the author has exploited 
certain properties of the operator 0 introduced in Section 4 (see equations (4-69.3) 
and (4-69.4) and the definition following). The properties alluded to are among those 
expressed by the relations 

0 (-A) = -0 (A) = [£(A)] T 

0 (A ± B) = 0 (A) ± 0 (B) 

0 (A) 0 (B) = [0 (B) 0 ( A) ] T 

N N 

$ (A) 53 a. 0 (Eb) =53 a. 0 (A) 0 (Bp 

i=l i=l 
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dt {$ [(A)( B )] } - $ 
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(B) | 


A dm = J £ (A) dm 


j & [A (r.,t)]} = $ 
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te) ! “-"I 


It is not difficult to show that 
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l l 


i ^,11 

5 (I V i} ^0 = (ll V (i) =n ^ ^ ^ ^ n j V^jK ’ (E_1) 


rp -0 


j=l 


j=l K=1 


i = 1, ... NA, 


use being made of equations (4-69.3) and (4-69.4). The notation in equation (E-l) 
should be self explanatory , that is , one should strongly suspect that the symbol 
^ ^->(i) simply denotes the inertia matrix of appendage i, referred to x^y.Zj, when 

the appendage is in its undeformed state which it assumes when n^=(f. The matrix 

(I 1 ) m will, for each i, be a known matrix of constants (provided by the manu- 

tt -0 

facturer of the appendage). 

Having found I 1 via equation (E-l), the matrix Q 1 is determined by 

Q 1 = I 1 - m. $ (tS 0 -* + n^ 1 ' ) ) & (-£^ o) - n^) , (E-2) 

[~| being the inertia matrix of appendage i referred to axes x!y!z! which have origin 
at the instantaneous CM of the appendage and are oriented as 

In writing equation (E-2), a direct application was made of what might be 
called the "generalized transfer theorem for inertia matrices" which is expressible as 

0 X yz =e ^ M £(q> £<-q> (e-3) 
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wherein [] ^ I s the inertia matrix of an arbitrarily shaped body of mass M referred 

to right-handed rectangular axes with origin at the CM of the body which has 
position vector q relative to the right-handed rectangular axes xyz; -^7 = f7(xyz-*- £n O 
denotes the rotation matrix defining the transformation from the xyz vector basis to 
the £ vector basis; and 0 X yz * s tlie * nert i a matrix of the body referred to the 
axes xyz. The expression £(q) $(-q) in equation (E-3) replaces its equivalent, 

(q T q) I( 3x3 ) - qq T , of the author's previous work (Reference 1), the symbol 
1(3x3) den °ting the 3x3 identity matrix. Admittedly, use of the expression is 
the more convenient. 

Repeated application of equation (E-3) leads to the following equation for the 
system inertia matrix [] (referred to the B-frame) 


NA 


i f + 7° 

CM 

+ m 0 £(%) £<-%) + Z D 1 T i + m i ^(qp 
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NSE 
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+ E 
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NR 

NSDOF 

+ e 

{T Ri lRiT Bi +m Ri*3Ri>S<-5 Ri )> + E C4 I® T Gi 

i=l 

i=l 

+ m Gi 

5(q' Gi ) «-q Gi ) + T P , I* Tgl + *3^) SC^)) 

N2DOF 

♦ E 

{T BOG 1 T BOG + m OG ^ ( ^OG ) ^ (_ ^OG ) + T BIG lIG T BIG 


i=l 


+ m^ ^(Qjq) &( Qjq) + T gg 1^ T jjg + m g ^^lg) q g^(i) * ( E 

By direct differentiation of equation (E-4), there follows, after further manipu- 
lation (giving due regard to the definitions and relations in Section 3 and to the 
relevant properties of the operator 0 ) , 


119 


[]= I f + m Q {£(q o ) £(-q o ) + jS(q Q ) ^(-q Q )} 


NA 


+ L <$l Q* 5 i + T, T (SfQ 1 -0‘ 5?) T. + m. [*(£,) Jgc-qj) + 4(5 j) £(-^>1) 
i=l 


NSE 

+ Z I Ei 


r Ei ~T 


i=l 


Ei (il Ei 1 " 1 fi Ei } T Ei + m Ei t ^ (q Ei ) ‘^ ( ' q Ei ) + ^ (q Ei ) ^ ( ~ q Ei )]} 


NP 


NR 


+ Z m pi »(q pi ) 4<-q pj > + A(q pl ) 4<-q p .)> + £ m Ri «(q„.) *(-q_.) 


i=l 


i pi ' i.. Ri I^V4 Ri > A>v M Ri > 

i=l 


NSDOF 

+ £( 5 Ri ) *(- 5 Ri ) } + £ ((A, i Gi 


i=l 


I G * ) T 
ui ' A Gi 


+ m Gi i^(q Gi ) ^(-q G i ) + ^ (q ci ) ^ ( ' q Gi )] + T ^i lgi - lgl q L) T , 


‘Gi'~"*Gi^ -gi '"gi gr-gi 


+ T Gi (n Gi lS1 ■ ,gI a Gi ) T Gi + m ii [ ^ < ?gi ) £<-3gi> + £<5 „h> *<-$„, 
N2DOF 


l gi M gi 


)]} 


Y {t t rn T n OG + t ig + t t T g T ^ 

^ U BOG l ^OG (I 1 T OGIG 1 T OGIG ) 

i=l 


(I GG + I^ G + tS t 'l o T IT + T^ rn^ /- T IG j. tST'* 

OGIG 1 T OGIG } Q OG ] T BOG T BIG [fJ I G (I I&) 


d IG + I*) af 0 j t b[g + t l «£ is - jg ah T 


B S g g J Bg 


m OG ^ ( 3 0 g> *-?og> + ^q 0G ) 


+ m iG ^Hq) A-?i G ) + ^(q IG ) *<-q IG » 


+ m g tf<q g > *<-? g > + S(-L»i„ 


g (i) 


(E-5) 


120 


Regarding the matrix i, all that will be said is that it is supposed here that a 
routine method exists for its computation.* As for the time derivatives of the position 
vectors , it is easily seen that they are given by 


q i = - r CM + ‘i 


T? + n? u/ 0 -* + 


i i 


q Ei " -r CM ■ £ Ei A Ei 


q Pi r CM + 5 Pi A Pi 


q 5 _r CM 


, 5 = 0, Ri, Gi, gi, OG, IG , g. 


while from equations (E-l) and (E-2), one has 
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N. N. 

l l 
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^ rj + ^rj ) ■ Sf HL n K + n j n K ) ^jK 


1=1 


j=l K=1 
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V = i 1 + m. n (l) ) £U. (o) + n (l) ) 


+ £(t (o) + ¥ (i) n (i) ) ^('F (i) ^ (i) )> • 


* Among the engineers engaged in digital simulation, few, if any, will be inclined to 
retain the term F in the expression for [] , not to mention certain of the other 
terms. In fact, most will modify the moment equation by deleting the term □ W B * 
One who is reluctant to retain all the terms belonging to [jwg should be cautioned 
not to discard terms which could be of major importance, in particular the CMG 
terms. 
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APPENDIX F 


ON TIME INDEPENDENT INTEGRALS 

Appearing in the moment equation, the bending equations, and the equation 
corresponding to the system coordinate 0. are terms having as a factor one of the 
time independent integrals defined below, their independence of time being a conse- 
quence of the assumption that the mode shape functions \ j = 1, ... N., are 

invariant under a rotation of flexible appendage i relative to the rigid central carrier. 
In a simulation program based upon the equations of this paper (and pertinent to a 
specific vehicle), those integrals would need to be evaluated but once. 
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